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1.  INTRODUCTION  AND  SUMMARY 
1,1  Introduction 

The  theory  of  Compressive  Sampling  has  the  potential  of  aehieving  a  reduetion  in 
the  eommunieation  and  eomputation  eosts  at  the  sensor.  A  discussion  of  the  theory  of 
compressive  sensing  is  in  section  2.1,  while  in  section  2.2  we  discussed  about  the 
research  work  carried  out  to  determine  whether  it  is  possible  to  predict  the  best  basis  for 
effective  reconstruction  of  different  types  of  signals  using  compressive  sensing 
algorithms.  The  performance  of  the  Nyquist  Folding  analog-to-information  receiver 
(NYFR),  a  compressive  sampling  based  application  for  signal  detection,  for  wideband 
signal  detection  in  multipath  fading  environment  is  discussed  in  section  2.3,  while 
another  compressive  sensing  based  method  for  spectrum  sensing,  which  makes  use  of 
signal  reconstruction  is  discussed  in  section  2.4. 

In  section  3,  upper  bounds  on  the  link  spectral  efficiency  of  amplify-and-forward 
cooperative  diversity  networks  with  independent  but  non-identically  distributed  wireless 
fading  statistics  are  studied  by  deriving  the  ergodic  capacity  of  three  distinct  adaptive 
source  transmission  techniques:  (i)  constant  power  with  optimal  rate  adaptation  (ORA); 
(ii)  optimal  joint  power  and  rate  adaptation  (OPRA);  and  (iii)  fixed  rate  with  truncated 
channel  inversion  (TCIFR).  Asymptotic  capacity  bound  is  also  derived  which  shows  that 
optimal  rate  adaptation  with  constant  power  policy  provides  roughly  the  same  ergodic 
capacity  as  the  optimal  joint  power  and  rate  adaptation  policy  at  high  mean  signal -to- 
noise  ratios  (SNRs).  Different  from  previous  related  studies,  we  advocate  a  simple 
numerical  procedure  for  unified  analysis  of  ergodic  channel  capacity  in  a  myriad  of 
fading  environments.  This  framework  allows  us  to  gain  insights  as  to  how  fade 
distributions  and  dissimilar  fading  statistics  across  the  diversity  paths  affect  the 
maximum  transmission  rates,  without  imposing  any  restrictions  on  the  fading  parameters. 
Analytical  expressions  are  also  derived  for  the  achievable  spectral  efficiency  and  average 
bit  error  rate  for  M-QAM  transmission  with  amplify-and-forward  cooperative  diversity 
network  in  a  Nakagami-m  fading  environment.  The  analysis  is  based  on  upper  bound, 
lower  bound  and  tight  approximation. 
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2,  COMPRESSIVE  SENSING  APPLICATION  IN  SPECTRUM  ACCESS 


2.1  Methods  Assumptions  and  Procedures  for  Compressive  Sensing  Applications 

One  of  the  reeent  researeh  areas  in  the  area  of  signal  proeessing  is  Compressive 
Sampling  or  Compressed  Sensing  as  it  may  be  referred  to  by  other.  This  theory  tends  to 
enable  a  reduetion  in  the  communication  and  computation  costs  at  the  sensor.  According 
to  Candes,  Tao  and  Romberg  [1],  a  small  number  of  random  projections  of  a  signal  that  is 
compressible  is  all  the  information  required  in  order  to  effectively  and  adequately 
reconstructs  the  signal.  In  other  words,  the  computational  burden  is  shifted  to  the  receiver 
side  of  the  communication  system  and  will  also  result  in  efficient  usage  of  the 
transmission  or  communication  bandwidth. 

Signal  compression,  in  the  traditional  sense  of  it,  heuristically  we  know  of  bases  or 
frame  that  can  be  used  to  model  a  particular  type  of  signal.  In  this  sense,  compression  for 
a  lot  of  signals  is  possible;  simply  by  specifying  the  largest  coefficients  in  the  transform 
domain  and  using  this  set  of  coefficients,  a  good  reconstruction  is  achievable. 
Compressed  sensing  on  the  other  hand,  illustrates  that  an  A-sample  signal  that  is  .^T-sparse 
requires  just  a  number  of  projections  of  the  order  0(K  log(N))  along  random  directions  in 
order  to  reconstruct  the  signal  in  good  quality  with  high  probability.  By  .^T-sparse  we 
meant  that  the  sparse  representation  of  the  signal  has  K  non-zero  elements.  In 
compressive  sampling  theory,  the  largest  coefficients  are  not  measured  directly.  Instead, 
the  projections  of  the  signal  onto  a  second  set  of  basis  functions  are  measured.  These  are 
referred  to  as  the  measurement  basis.  The  two  bases,  sparsifying  basis  and  the 
measurement  basis  are  required  to  be  incoherent.  One  of  the  fundamental  theorems  of 
compressive  sampling  is  stated  below  [1]; 

Theorem  1:  Assume  that  a  signal  x  is  K-sparse  and  that  we  are  given  P  Fourier 
coejficients  with  frequencies  selected  uniformly  at  random.  Suppose  that  the  number  of 
observations  obeys 

P>CK\ogN.  (1) 

Then  minimizing  li  reconstructs  x  exactly  with  overwhelming  probability.  In  details,  if  the 
constant  C  is  of  the  form  22  (d  +  1)  in  Equation  (1),  then  the  probability  of  success 
exceeds 

1-  0(N-d). 

Considering  the  above  stated  theorem,  we  can  conclude  that  it  is  possible  to 
reconstruct  signal  x  without  suffering  any  information  loss,  just  by  taking  any  set  of  P 
frequency  coefficients.  Also  we  can  recover  signal  x  by  minimizing  a  convex  functional 
which  does  not  assume  any  knowledge  about  the  number  of  nonzero  coordinates  of  x, 
their  locations,  and  their  amplitudes  which  we  assume  are  all  completely  unknown  a 
priori.  According  to  Candes  [5],  P  is  in  general  the  optimal  number  of  coefficients 
required.  There  are  situations  whereby  the  minimum  number  of  samples  required  for 
exact  reconstruction  regardless  of  the  method  must  be  about  K  log  N . 

2.1.1  Reconstruction  (Sparsifying)  Bases 

Transform  bases  can  be  used  to  make  the  signal  sparse.  As  mentioned  earlier,  the 
reconstruction  bases  that  are  considered  for  this  work  includes  Fast  Fourier  Transform 
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(FFT),  Discrete  Cosine  Transform  (DCT)  and  Diserete  Wavelet  Transform  (DWT)  with 
Haar,  Daubeehies  4,  Coiflet  and  Symlet  filters. 


Fast  Fourier  Transform 

The  Fourier  transform  is  widely  used  over  linear  time-invariant  signal  proeessing 
beeause  sinusoidal  waves  e“"‘are  eigenvectors  of  linear  time-invariant  operators.  It 
involves  the  transformation  of  signal  waveforms  into  an  infinite  summation  (the  integral) 
of  sinusoids  with  the  frequency  differences  between  adjacent  sinusoids  infinitesimally 

small.  Fourier  Transform  ean  be  expressed  as  [6] 

00 

5(/)- (2) 


where  S(f)  is  the  Fourier  transform  of  the  signal  s(t),  f  is  the  frequeney  and  t  is  time. 
The  Fourier  transform  provides  simple  answers  to  most  questions.  For  a  signal  s(t)  of  N 
points,  a  direct  calculation  of  the  N  diserete  Fourier  sums  is  defined  by; 

-ilTrkn'^ 


N-l 


n=0 


N 


for  0  <k<N , 


(3) 


The  Fast  Fourier  Transform  (FFT)  algorithm  reduees  the  numerieal  complexity  to 
0(N  log2  N)  by  reorganizing  the  calculations.  Fast  Fourier  Transform  (FFT)  has  been 
used  extensively  for  signal  analysis  in  the  past  and  is  still  being  used.  Frequency 
resolution,  phase  measurements  and  error  calculations  are  some  of  the  applieation  areas 
of  FFT.  It  has  the  ability  to  localize  signals  in  the  frequeney  domain. 


Discrete  Cosine  Transform 

Discrete  Cosine  Transfroms  (DCT)  has  been  applied  in  the  past  in  the  area  of  data 
compressions,  feature  extraetion,  multiframe  deteetion,  etc.  Basis  vectors  in  the  DCT 
domain  have  been  proved  to  provide  good  energy  compaction  for  natural  images.  Given 
an  N-dimensional  vector  input  x,  the  N  point  DCT  Y,  is  defined  as 

Y  =  D^.x  (4) 


where  Dn  is  the  DCT  transformation  matrix  of  size  N  by  N,  with  elements  DN(i,j)  defined 
as; 


c, 


(2j  +  V)i7r 
2N 


(5) 


where 


|l/>/2  i-O 
|l  i>0 


The  DCT  is  orthogonal,  hence  the  inverse  transformation  ean  be  written  as 


x  =  D‘^.Y 


(V) 


Wavelet  Transform 

Wavelet  transform  has  the  capability  of  localizing  signals  in  both  time  and  frequency 
domains.  As  a  result  of  this  speeial  eharacteristie,  it  has  found  usage  in  various 
applieations  such  as  image  proeessing,  biomedieal  image  system,  electromagnetic 
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systems,  data  compression,  signal  detection  [7,  8],  and  so  on.  By  definition,  a  wavelet  \|; 

is  a  function  of  zero  average  [6]: 

00 

^y/{t)dt-0  (8) 


which  is  dilated  with  a  scale  parameter  5,  and  translated  by  u: 


1 

ds 


''  t -u'^ 


(9) 


where  y/it)  is  the  mother  wavelet  equation,  u  is  the  translation  parameter  and  s  is  the 
dilation  factor  also  referred  to  as  the  scale.  The  wavelet  transform  of  a  function  or  signal 
/,  at  the  scale  5  and  position  u  is  computed  by  correlating /with  a  wavelet  atom  as  shown 
in  equation  9. 


OO  . 

Xs 

where  Wf(u,s)  is  the  wavelet  transform,  f(t)  is  the  original  signal  [6].  The  integral 
function  represents  the  Continuous  Wavelet  Transform  (CWT). 


t -u 


dt 


(10) 


2.1.2  Measurement  (Sampling)  Bases 

One  of  the  goals  of  compressive  sampling  signal  reconstruction  is  to  be  able  to  design 
a  matrix  O,  which  satisfies  the  requirement  of  orthogonality  for  any  subset  of  columns  of 
size  K.  Though  it  might  be  difficult  to  prove  that  there  exists  a  matrix  that  will  obey  UUP 
for  large  values  of  K,  we  know  that  trivial  randomized  construction  will  do  so  with  good 
probability.  Examples  of  such  random  matrices  that  are  capable  of  achieving  results  are 
discussed  below: 

a)  Gaussian  measurements:  The  sensing  matrix  O  is  obtained  by  taking  independent 
samples  from  the  normal  distribution  with  zero  mean  and  a  variance  of  value  1/P. 
Then  if 

K<C.  P/log(N/P)  (11) 

[9],  [10]. 

b)  Binary  measurements:  In  this  case,  the  sensing  matrix  O  is  obtained  by  taking 
independent  samples  from  the  symmetric  Bernoulli  distribution  with  parameter; 

Prob(0  =  ±1  /  ^/P)  =  1  /  2  (12) 

c)  Fourier  measurements:  Let  us  assume  that  the  sensing  matrix  O  is  a  partial  Fourier 
matrix  obtained  by  taking  P  rows  uniformly  at  random,  but  renormalizing  the 
columns  such  that  they  are  unit-normed.  Candes  and  Tao  [11]  have  shown  that  CS 
theorem  holds  with  overwhelming  probability  if 

K<C.  P/(logA0'^  (13) 

This  result  was  recently  improved  by  Rudelson  and  Vershynin  [12]  to 

K<C.  P/(logN)^  (14) 

d)  Incoherent  measurements:  Let  us  assume  now  that  the  sensing  matrix  O  is  obtained 
by  selecting  P  rows  uniformly  at  random  from  an  N  by  N  orthonormal  matrix  U  and 
renormalizing  the  columns  so  that  they  are  unit-normed.  U  can  be  considered  as  the 
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matrix  O^F*  which  maps  the  object  from  the  'F  to  the  O-domain.  The  argument  used 
to  prove  that  UUP  holds  for  the  Fourier  measurements  can  be  extended  to  this  case  as 
well.  CS  theorem  holds  with  high  probability  provided  that  Equation  (15)  is  satisfied. 


K<C~ 


P 

(\ogNY 


(15) 


where  u  :=  ^//V  max.  U,  ■ 


2.1.3  Minimization  Algorithms 

Minimization  algorithm  methods  for  CS  come  in  two  different  categories  namely; 

i.  Greedy  methods  -  makes  a  sequence  loeally  optimal  choice  in  an  effort  to  determine 
a  globally  optimal  solution.  Example  includes  Matching  Pursuit  and  Orthogonal 
Matching  Pursuit. 

ii.  Convex  relaxation  methods  -  replaces  a  combinatorial  sparse  approximation  problem 
with  a  related  convex  program.  Example  includes  Basis  Pursuit.  The  minimization 
algorithms  are  discussed  below: 


Matching  Pursuit:  If  the  dictionary  is  orthonormal,  the  sparse  approximation  problem 
admits  a  straightforward  algorithm.  It  is  possible  to  build  a  solution  one  term  at  a  time  by 
selecting  at  each  step  the  atom  that  eorrelates  most  strongly  with  the  residual  signal. 
Matching  pursuit  (MP)  extends  this  idea  to  other  types  of  dietionary.  MP  begins  by 
setting  the  initial  residual  equal  to  the  input  signal  and  making  a  trivial  initial 
approximation.  That  is, 

ro=s  and  ao  =  0. 

At  step  k,  MP  ehooses  another  index  kk  by  solving  an  easy  optimization  problem 


4  Gargmax(r,_j,(j9^) 

Then  it  calculates  a  new  approximation  and  a  new  residual 

(16) 

~  %-i  "*■  (^-1  ’  ^^4  )  ^^4 

and 

(17) 

The  residual  can  also  be  expressed  as 

(18) 

r,=s-a, 

(19) 

When  the  dictionary  is  an  orthonormal  basis,  the  approximant  is  always  an  optimal-term 
representation  of  the  signal.  MP  is  considered  as  the  pure  greedy  algorithm  [13]. 


Orthogonal  Matching  Pursuit:  OMP  adds  a  least-squares  minimization  to  eaeh  step  of 
MP  to  obtain  the  best  approximation  over  the  atoms  that  have  already  been  ehosen.  This 
revision  significantly  improves  the  behavior  of  the  algorithm.  OMP  is  initialized  the  same 
way  as  MP,  and  at  eaeh  step,  a  new  atom  is  selected  according  to  the  same  rule  as  MP, 

via  Equation  (17).  But  the  approximants  are  calculated  differently.  Eet  Ak  =  [ki . kj 

list  the  atom  that  has  been  ehosen  at  step  k.  Then  the  kth  approximant  is 
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l2 


(20) 


def 

a^_  —  argminlU  — a 

a 

subject  to  a  e  span  :  A  ^  } 

This  minimization  can  be  performed  inerementally  with  standard  least-squares 
teehniques.  As  before,  the  residual  is  ealculated  as, 

rk=s-a,. 

Note  that  OMP  never  selects  the  same  atom  twiee  beeause  the  residual  is  orthogonal  to 
the  atoms  that  have  already  been  ehosen.  Consequently,  the  residual  must  equal  zero  after 
d  steps. 

Basis  Pursuit:  Convex  relaxation  offers  another  approaeh  to  sparse  approximation.  The 
fundamental  idea  is  that  the  number  of  terms  in  a  representation  ean  be  approximated  by 
the  absolute  sum  of  the  eoeffieients.  This  absolute  sum  is  a  eonvex  funetion,  and  so  it  ean 
be  minimized  in  polynomial  time.  Given  an  input  signal  s,  the  BP  problem  is 

minV  1^  I  (21) 

subject  to  \b^(p^  =  s\  (22) 

where  b^  is  a  eollection  of  eomplex  eoeffieients.  One  hopes  that  the  nonzero  eoeffieients 
in  the  solution  of  the  BP  problem  will  identify  the  atoms  in  the  optimal  representation  of 
the  input  signal  and  their  eoeffieients. 

Least  Absolute  Shrinkage  and  Selection  Operator  (LASSO) 

LASSO  (Least  Absolute  Shrinkage  and  Selection  Operator)  is  derivatives  of  LARS 
(Least  Angle  Regression)  minimum  fi  norm  reeonstruction  algorithm.  This  greedy 
algorithm  take  advantage  of  geometrie  eorrelation  between  the  targeted  and  the 
eoeffieients  most  correlated  with  the  measured  response.  LASSO  algorithms  take 
advantage  of  quadratie  programming  eoneepts  beeause  they  produce  a  weighted 
minimum  f  i  and  £2  norm  solution.  The  LASSO  algorithm  speoifieally  solves 

mini  ||F-5||  -l-/l||5||i  sueh  that  T  =  05  (23) 

[2m II  II2  II  IIJ 

for  any  parameter  X.  The  value  of  X  typieally  eorresponds  to  the  maximum  eorrelation 
between  the  measurements  Y  and  the  observation  matrix  O. 

Stepwise  Orthogonal  Matching  Pursuit  (StOMP) 

For  n  large.  Stepwise  Orthogonal  Matehing  Pursuit  (StOMP)  correetly  reeovers  the 
sparsest  solution  of  Y  =  Oy  over  a  region  of  the  sparsity/indeterminaey  plane  eomparable 

to  the  region  where  \i  minimization  is  suceessful.  StOMP  sueeessively  transforms  the 
signal  into  a  negligible  residual.  Starting  with  initial  residual  ro  =  Y,  at  the  s-th  stage  it 
forms  the  ‘matehed  filter’  O^rs-i,  identifies  all  eoordinates  with  amplitudes  exceeding  a 
speeially-ehosen  threshold,  solves  a  least-squares  problem  using  the  seleeted  eoordinates, 
and  subtraets  the  least  squares  fit,  produeing  a  new  residual.  After  a  fixed  number  of 
stages  it  stops.  In  eontrast  to  Orthogonal  Matehing  Pursuit  (OMP),  many  eoeffieients  ean 
enter  the  model  at  each  stage  in  StOMP  while  only  one  enters  per  stage  in  OMP;  and 
StOMP  takes  a  fixed  number  of  stages  ,  while  OMP  can  take  many  StOMP  runs  mueh 
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faster  than  competing  proposals  for  sparse  solutions,  such  as  '  1  minimization  and  OMP, 
and  so  is  attractive  for  solving  large-scale  problems  [14], 

2.2  Sparseness  Measure  of  Signals  for  Compressive  Sampling 

Choice  of  reconstruction  or  sparsity  basis  is  one  of  the  major  criteria  in  order  to  achieve  a 
good  CS  signal  reconstruction.  One  of  the  goals  of  this  research  work  is  to  be  able  to 
show  how  different  choice  of  sparsity  basis  affects  the  quality  of  the  reconstructed  signal. 
One  of  the  objectives  of  this  work  is  to  be  able  to  adequately  predict  the  best  basis. 

2.2.1  Sparseness  Measures 

Sparse  signal  representation  helps  a  great  deal  in  achieving  robust  methods  for 
compression,  detection,  denoising,  and  signal  separation  [15,  16].  The  concept  of 
sparsity  has  been  used  to  find  meaningful  data  representation  in  the  areas  of  image 
analysis,  vision  research,  speech  and  music  analysis  [17,  18,  19].  In  simple  terms, 
sparsity  of  a  signal  can  be  defined  as  the  number  of  non-zero  elements  in  a  vector 
representation  of  the  signal.  This  implies  that  a  signal  that  is  considered  sparse  has  most 
of  its  components  to  be  zero.  In  a  practical  sense  however,  sparsity  means  that  most 
signal  components  are  small  enough  not  to  contribute  significantly  to  the  information 
content  of  the  signal. 

Some  of  the  sparseness  measures  found  in  literature  [16]  [16,  20]  are  described  in  Table 
1  and  discussed  below.  A  common  measure  is  the  f  pseudo-norm  given  by  the  number 
of  non-zero  elements  as  described  earlier.  The  f  measure  does  not  take  into 
consideration  the  magnitude  of  the  non-zero  element.  Hence  replacing  a  non-zero 
element  with  another  one,  does  not  affect  the  sparseness  measure.  This  does  not  give  a 
true  representation  of  the  distribution  of  the  component.  In  the  presence  of  noise,  f(s)  is 
used  because  the  noise  results  in  very  few  of  the  components  being  exactly  zero. 


Table  1  Sparseness  Measures. 


l\s) 

' — 

IV 

/' 

1” 

M 

■k 

O 

A 

Log 

Xdog(l+x^.f 

tanh(a,b) 

Z/anh(flx/ 

Kurtosis 

(k4) 

y  ixJ' 

^j\  d 

(U,\1 

Gini  Index 

Discussed  in  text 
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Adding  noise  to  a  signal  can  change  the  sparseness  measure  when  using  f  norm 
measure,  hence  f(e)  is  used  instead.  The  F  norm  is  often  used  instead  of  f(e)  because  it 
has  a  useable  gradient  and  thus  NP-hard  combinatoric  optimization  problems  can  be 
avoided.  The  tanh(a,b)  measure  is  sometimes  used  instead  of  F,  p<l  norm  because  each 
component  contribution  is  limited  to  range  (0,  1).  The  concave  property  of  both  of  them 
enforces  sparsity.  The  log  measure  is  concave  outside  some  range  but  convex  near  the 
origin.  In  all  of  the  sparseness  measures  described  in  Table  1  with  the  exception  of 
kurtosis  and  the  Gini  index,  lower  values  correspond  to  more  sparsity. 

The  Gini  index  is  a  measure  used  by  economists  in  evaluating  the  inequality  of  wealth 
distribution.  Dalton  defines  four  properties  that  a  sparse  measure  should  satisfy  [21]  and 
the  Gini  index  satisfies  all  conditions  [20].  To  define  the  Gini  index  for  a  real-valued 
data  set  {Z},  let  x  =  {xi,  X2....xn}  represent  the  data  elements  of  {Z}  ordered  in  ascending 
absolute  value  such  that  |xj<|;c2|<|;c3|<  -  |x^|,  and  define  a  piecewise  linear  function  with 

values  between  0  and  1  as 

=  fori=l,....iV  (24) 

where  L(  )  is  denoted  as  the  Lorenz  curve.  The  Lorenz  curve  equation  can  be  rewritten  as, 

Shi 

^  for/=l,....A  (25) 

where  is  the  F  norm  sparseness  measure. 

The  Lorenz  curve  can  be  likened  to  the  cumulative  distribution  plot  of  a  statistical 
data;  it  does  give  a  good  representation  of  how  sparsely  distributed  the  data  is.  The  Gini 
index  is  twice  the  area  between  the  Lorenz  curve  and  the  45  degree  line.  The  area 
underneath  the  Lorenz  curve  is 

=  +  (26) 

^2''  n=\ 

Therefore,  the  Gini  index  is  defined  as 

G(x)  =  l-2A(x)  (27) 

In  economic  parlance,  the  distribution  where  all  individuals  have  equal  wealth  is  the 
least  sparse,  whereas  the  distribution  where  one  individual  controls/possesses  all  the 
wealth  is  the  sparsest.  In  a  similar  manner,  the  data  distribution  whereby  all  the  data  have 
equal  value  is  the  least  sparse,  whereas  the  one  where  there  is  one  big  value  and  the  rest 
are  zero,  is  the  sparsest. 

The  value  of  the  Gini  index  G{x)  always  lies  between  0  and  1 ,  with  1  being  the  sparsest 
and  0  representing  no  sparsity.  One  of  the  advantages  of  this  sparseness  measure  is  that 
normalization  is  not  required  when  comparing  vectors  of  different  lengths;  for  the  other 
sparseness  measures  discussed  in  this  section,  some  form  of  normalization  is  required 
before  comparing  vectors  of  different  lengths.  Also  the  Gini  Index  gives  an  adequate 
sparseness  measurement  in  the  presence  of  noise  unlike  the  f  norm  which  treats  every 
value  other  than  zero  as  the  same. 
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2.2.2  Results  and  Discussion  of  1-D  Compressive  Sampling  Signal  Reconstruction 
without  Noise 

Six  different  types  of  signals  were  simulated,  namely  Heavisine,  Bump,  Blocks,  Ramp, 
Doppler  and  Linear  chirp.  The  bases  considered  are  Fast  Fourier  Transform  (FFT), 
Discrete  Cosine  Transform  (DCT),  and  different  families  of  Discrete  Wavelet  Transform 
(DWT).  We  investigated  how  sparse  each  of  the  signals  is  while  represented  in  the  bases 
being  considered  using  the  various  sparseness  measures  listed  in  Table  1.  Thereafter,  CS 
reconstruction  was  carried  out  to  determine  if  there  is  a  relationship  between  the  best 
basis  choice  of  the  sparseness  measures  and  the  mean  square  error  of  the  CS 
reconstruction  carried  out. 

The  signals  were  considered  with  varying  signal  lengths  of  128,  256,  512,  1024,  2048 
and  4096.  Figure(s)  1  to  7  show  the  plot  of  the  sparseness  measures  for  the  block  signal 
type,  at  different  original  signal  lengths.  It  can  be  seen  from  the  plots  that  for  the  block 
signal  type,  all  the  sparseness  measures  except  kurtosis,  predicted  that  the  best 
sparsifying  basis  is  the  haar  wavelet. 


Lg  norm  of  Block  Signal  vs  Signal  Length 


Figure  1  LO  norm  measures. 
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Norm 


Lp  norm,  p  =  0.5  of  Block  Signal  vs  Signal  Length 


Figure  2  Lp  norm  measure,  p  =  0.5. 


L.|  norm  of  Block  Signal  vs  Signal  Length 
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Log  Measure 


Log  Measure  of  Block  Signal  vs  Signal  Length 


Tanh  of  Block  Signal  vs  Signal  Length 
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Kurtosis  of  Block  Signal  vs  Signal  Length 


Gini  Index  of  Block  Signal  vs  Signal  Length 


Similar  data  were  obtained  for  the  other  type  of  signal  that  was  considered.  Tables  2 
to  7,  contains  the  data  obtained  for  all  the  signal  type  considered.  Tables  2  to  5  illustrated 
the  sparsifying  basis  predicted  by  each  of  the  sparseness  measures  for  each  of  the  signal 
type  considered  at  different  signal  length.  However,  two  particular  cases  of  signal  lengths 
2048  and  4096  were  considered  and  the  results  obtained  were  tabulated  in  Tables  6  to  7. 
For  these  two  signal  lengths,  the  sparseness  of  the  signals  was  measured  and  the  signals 
were  reconstructed  using  compressive  sampling  at  a  decimation  ratio  of  4:1.  The 
reconstruction  was  carried  out  using  basis  pursuit  and  orthogonal  matching  pursuit.  Mean 
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square  error  of  the  reeonstrueted  signals  in  eomparison  to  the  original  signals  were 
eomputed  and  used  as  performanee  metrie.  Apart  from  Kurtosis  which  was  mostly  wrong 
in  its  prediction,  the  other  measures  performed  averagely  well  with  Gini  index  and  lo 
norm  having  the  best  prediction  across  board. 


Table  2  Sparseness  Measure  of  Various  1-D  signals  (signal  length  =  128). 


FFT 

DCT 

WAV 

ELET 

Haar 

Coiflet 

Db4 

Symiet 

Lq 

0.6172 

0.0547 

0.1797 

0.0859 

0.0859 

0.0547 

Lp(p=l) 

8.9261 

0.8684 

1.0995 

0.8688 

0.9144 

0.8630 

Lp(p=0.5) 

482.35 

43.64 

77.712 

18.803 

22.684 

18.402 

Log 

0.9157 

0.1577 

0.2423 

0.1625 

0.1833 

0.1628 

Tanh 

0.9157 

0.2471 

0.4423 

0.1633 

0.1933 

0.1618 

Kurtosis 

0.3511 

0.2787 

0.3812 

0.1398 

0.1577 

0.1539 

Heavisine 

Gini  Index 

0.8017 

0.8300 

0.7261 

0.9176 

0.9083 

0.9188 

Lq 

0.2500 

0.1094 

0.1172 

0.1328 

0.1484 

Lp  (p=l) 

0.5714 

0.3607 

0.3807 

0.4016 

0.3901 

Lp(p=0.5) 

828.48 

61.291 

19.042 

25.097 

26.648 

24.484 

Log 

1.5714 

0.1450 

0.1049 

0.1010 

0.1098 

0.1095 

Tanh 

0.9848 

0.5729 

0.2535 

0.3012 

0.3105 

0.2918 

Kurtosis 

0.0317 

0.0376 

0.0704 

0.0912 

0.0657 

0.0708 

Bumps 

Gini  Index 

0.3142 

0.4224 

0.7694 

0.7098 

0.7005 

0.7204 

Lq 

0.3047 

0.1250 

0.2266 

0.1875 

0.1953 

Lp(p=l) 

1.2464 

0.9151 

1.1021 

1.1401 

1.1173 

Lp(p=0.5) 

1468.2 

108.93 

25.253 

68.707 

65.362 

70.559 

Log 

1.9230 

0.3296 

0.2452 

0.2959 

0.3111 

0.3037 

Tanh 

0.9910 

0.6616 

0.2399 

0.4515 

0.4471 

0.4569 

Kurtosis 

0.1889 

0.2058 

0.1918 

0.2053 

0.1848 

0.2080 

Blocks 

Gini  Index 

0.5465 

0.6087 

0.8539 

0.7337 

0.7435 

0.7202 

Lo 

0.6172 

0.0469 

0.0391 

0.0469 

0.0391 

0.0469 

Lp  (p=l) 

2.2908 

0.1806 

0.1604 

0.1097 

0.1274 

0.1121 

Lp(p=0.5) 

244.01 

17.670 

13.028 

4.1286 

6.2330 

4.5294 

Log 

0.6384 

0.0284 

0.0275 

0.0240 

0.0249 

0.0234 

Tanh 

0.9950 

0.1231 

0.1285 

0.1006 

0.1053 

0.0973 

Kurtosis 

0.0531 

0.0822 

0.1030 

0.1863 

0.1538 

0.2097 

Doppler 

Gini  Index 

0.4565 

0.5342 

0.6535 

0.8484 

0.7879 

0.8331 

Lo 

0.4531 

0.0313 

0.0078 

0.0156 

0.0156 

0.0234 

Lp(p=l) 

1.7133 

0.1445 

0.0958 

0.0859 

0.0841 

0.0852 

Lp(p=0.5) 

157.4500 

12.0380 

4.6317 

1.6738 

1.4626 

1.5895 

Log 

0.4041 

0.0264 

0.0228 

0.0225 

0.0220 

0.0224 

Tanh 

0.9316 

0.0858 

0.0633 

0.0638 

0.0692 

0.0646 

Kurtosis 

0.1668 

0.1527 

0.2350 

0.2411 

0.2774 

0.2613 

Ramp 

Gini  Index 

0.5686 

0.6417 

0.8413 

0.9223 

0.9254 

0.9235 

Lo 

0.7344 

0.3203 

0.2578 

0.1875 

0.2578 

0.1797 

Lp  (p=l) 

7.7159 

0.6115 

0.5353 

0.4285 

0.4526 

0.4519 

Lp(p=0.5) 

985.50 

71.0530 

55.8450 

32.5980 

37.6560 

36.9180 

Log 

1.7783 

0.1543 

0.1376 

0.1198 

0.1225 

0.1242 

Tanh 

1.0000 

0.7077 

0.5533 

0.3631 

0.4204 

0.3977 

Kurtosis 

0.0080 

0.0122 

0.0224 

0.0330 

0.0345 

0.0301 

Chirp 

Gini  Index 

0.0309 

0.2872 

0.4434 

0.6389 

0.5938 

0.6012 

Approved  for  Public  Release;  Distribution  Unlimited. 
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Table  3  Sparseness  Measure  of  Various  1-D  signals  (signal  length  =  256). 


FFT 

DCT 

WAV 

ELET 

Haar 

Coiflet 

Db4 

Symiet 

Lq 

0.6016 

0.0313 

0.I0I6 

0.0391 

0.0547 

0.0430 

Lp  (p=l) 

9.4953 

0.6508 

0.8223 

0.5928 

0.6583 

0.6173 

Lp  (p=0.5) 

923.39 

56.318 

93.769 

14.007 

18.821 

14.920 

Log 

0.8948 

0.1047 

0.I65I 

0.0950 

0.II49 

0.0991 

Tanh 

0.9202 

0.1593 

0.2706 

0.0954 

0.I2I7 

0.1026 

Kurtosis 

0.3515 

0.2779 

0.3796 

0.I7I6 

0.1585 

0.1442 

Heavisine 

Gini  Index 

0.8212 

0.8575 

0.8001 

0.9594 

0.9496 

0.9564 

0.7188 

0.1797 

0.0938 

0.I0I6 

0.II33 

0.0898 

Lp  (P=l) 

8.2609 

0.4740 

0.2979 

0.3079 

0.3099 

0.3158 

Lp  (p=0.5) 

1815.60 

96.201 

27.146 

30.532 

29.971 

31.632 

Log 

1.5997 

0.II36 

0.0844 

0.0845 

0.0866 

0.0879 

Tanh 

1. 0000 

0.4359 

0.2173 

0.2271 

0.2193 

0.2288 

Kurtosis 

0.0342 

0.0399 

0.0562 

0.0581 

0.0556 

0.0557 

Bumps 

Gini  Index 

0.4213 

0.4947 

0.8II3 

0.7901 

0.7931 

0.7844 

Lo 

0.7109 

0.2383 

0.0820 

0.1523 

0.1484 

0.1328 

Lp  (P=l) 

17.7490 

1.0280 

0.6957 

0.8575 

0.8791 

0.8498 

Lp  (p=0.5) 

3129.60 

164.95 

23.732 

74.103 

69.725 

73.673 

Log 

1.9342 

0.2512 

0.1720 

0.2178 

0.2216 

0.2152 

Tanh 

0.9996 

0.5365 

0.1539 

0.3156 

0.3199 

0.3220 

Kurtosis 

0.1878 

0.2043 

0.1905 

0.2016 

0.1827 

0.2084 

Blocks 

Gini  Index 

0.6098 

0.6639 

0.9120 

0.8254 

0.8332 

0.8219 

Lo 

0.5469 

0.0273 

0.0313 

0.0234 

0.0195 

Lp  (P=l) 

2.5805 

0.1438 

0.1278 

0.0837 

0.0940 

0.0807 

Lp  (p=0.5) 

488.93 

24.958 

18.197 

4.2194 

5.9734 

3.9437 

Log 

0.6222 

0.0240 

0.0200 

0.0204 

0.0193 

Tanh 

0.9705 

0.0859 

0.0624 

0.0757 

0.0635 

Kurtosis 

0.0533 

0.0850 

0.1409 

0.1545 

0.I7II 

Doppler 

Gini  Index 

0.5634 

0.6207 

■nBBW 

0.9041 

0.8741 

0.9090 

Lo 

0.4219 

0.0156 

0.0195 

0.0078 

0.0195 

Lp  (P=l) 

1.9359 

0.II47 

0.0692 

0.0620 

0.0603 

0.0627 

Lp  (p=0.5) 

326.96 

17.403 

4.9767 

1.3670 

1. 1357 

1.3434 

Log 

0.4III 

0.0210 

0.0174 

0.0168 

0.0168 

0.0173 

Tanh 

0.9316 

0.0564 

0.0400 

0.0436 

0.0401 

0.0429 

Kurtosis 

0.1667 

0.1527 

0.2343 

0.2401 

0.2770 

0.2642 

Ramp 

Gini  Index 

0.6174 

0.6839 

0.8924 

0.9571 

0.9608 

0.9572 

Lo 

0.7383 

0.3242 

0.2266 

0.1836 

0.2305 

0.1953 

Lp  (P=l) 

11.0300 

0.6189 

0.5344 

0.4286 

0.4387 

0.4386 

Lp  (p=0.5) 

2819.10 

144.25 

II0.I900 

63.3050 

70.4190 

66.3740 

Log 

2.0859 

0.1570 

0.1380 

0.II94 

0.II83 

0.1224 

Tanh 

1. 0000 

0.7140 

0.5357 

0.3739 

0.3953 

0.3872 

Linear 

Kurtosis 

0.0040 

0.0060 

0.0II6 

0.0189 

0.0267 

0.0179 

Chirp 

Gini  Index 

0.0241 

0.2832 

0.4550 

0.6438 

0.6I2I 

0.6290 

Approved  for  Public  Release;  Distribution  Unlimited. 
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Table  4  Sparseness  Measure  of  Various  1-D  signals  (signal  length  =  512). 


FFT 

DCT 

WAV 

ELET 

Haar 

Coiflet 

Db4 

Symiet 

Lq 

0.5918 

0.0273 

0.0664 

0.0254 

0.0293 

0.0215 

Lp  (p=l) 

10.0470 

0.4836 

0.6087 

0.4103 

0.4684 

0.4346 

Lp  (p=0.5) 

1793. 1 

74.011 

II0.I6 

I0.8II 

14.619 

11.437 

Log 

0.8837 

0.0686 

0.II25 

0.0575 

0.0708 

0.0601 

Tanh 

0.9192 

0.1027 

0.1769 

0.0587 

0.0712 

0.0598 

Kurtosis 

0.3520 

0.2769 

0.3775 

0.I87I 

0.1584 

0.1529 

Heavisine 

Gini  Index 

0.8348 

0.8734 

0.8545 

0.9786 

0.9736 

0.9775 

0.7090 

0.II52 

0.0762 

0.0703 

0.0820 

0.0723 

Lp  (P=l) 

10.1470 

0.4090 

0.2651 

0.2593 

0.2643 

0.2603 

Lp  (p=0.5) 

3996.4 

147.35 

39.5760 

35.775 

36.647 

36.486 

Log 

1.6058 

0.0983 

0.0763 

0.0734 

0.0752 

0.0730 

Tanh 

0.9970 

0.3198 

0.1675 

0.1683 

0.I7I8 

0.I7I7 

Kurtosis 

0.0330 

0.0389 

0.0459 

0.0503 

0.0480 

0.0510 

Bumps 

Gini  Index 

0.5356 

0.5950 

0.8543 

0.8574 

0.8547 

0.8542 

Lo 

0.7168 

0.1660 

0.0566 

0.1035 

0.0996 

0.0996 

Lp  (P=l) 

20.2330 

0.8229 

0.5204 

0.6463 

0.6579 

0.6413 

Lp  (p=0.5) 

6553.0 

239.22 

21.882 

74.505 

67.6290 

73.1190 

Log 

1.9430 

0.1846 

0.II60 

0.1543 

0.1535 

0.1537 

Tanh 

0.9997 

0.4022 

0.I0I5 

0.2I2I 

0.2163 

0.2II9 

Kurtosis 

0.I86I 

0.2035 

0.1878 

0.1985 

0.1804 

0.2073 

Blocks 

Gini  Index 

0.6570 

0.7098 

0.9485 

0.8890 

0.8967 

0.8880 

Lo 

0.4941 

0.0215 

0.0215 

0.0176 

0.0215 

0.0156 

Lp  (P=l) 

2.8054 

0.II09 

0.1002 

0.0623 

0.0687 

0.0604 

Lp  (p=0.5) 

895.96 

32.532 

24.326 

3.7409 

5.4810 

3.6227 

Log 

0.5707 

0.0196 

0.0187 

0.0160 

0.0163 

0.0155 

Tanh 

0.8450 

0.0622 

0.0612 

0.0408 

0.0543 

0.0409 

Kurtosis 

0.0537 

0.0843 

0.I06I 

0.II89 

0.1552 

0.1450 

Doppler 

Gini  Index 

0.6735 

0.7105 

0.7819 

0.9456 

0.9278 

0.9469 

Lo 

0.4434 

0.0II7 

0.0078 

0.0078 

0.0078 

0.0137 

Lp  (P=l) 

2.1556 

0.0900 

0.0503 

0.0455 

0.0435 

0.0462 

Lp  (p=0.5) 

670.98 

25.06 

5.4547 

1. 1547 

0.9322 

1. 1349 

Log 

0.4145 

0.0160 

0.0124 

0.0125 

0.0123 

0.0130 

Tanh 

0.9316 

0.0374 

0.0260 

0.0264 

0.0238 

0.0276 

Kurtosis 

0.1667 

0.1507 

0.2472 

0.2362 

0.2770 

0.2623 

Ramp 

Gini  Index 

0.6560 

0.7155 

0.9247 

0.9760 

0.9787 

0.9758 

Lo 

0.7383 

0.3203 

0.2266 

0.I7I9 

0.2070 

0.1758 

Lp  (P=l) 

15.7160 

0.6239 

0.5294 

0.4232 

0.4441 

0.4338 

Lp  (p=0.5) 

8037.3 

291.04 

215.33 

119.29 

137.51 

126.72 

Log 

2.3924 

0.1589 

0.I37I 

0.II93 

0.1226 

0.I2I6 

Tanh 

1. 0000 

0.7189 

0.5252 

0.3617 

0.3985 

0.3796 

Linear 

Kurtosis 

0.0020 

0.0030 

0.0068 

O.OIII 

0.0103 

0.0103 

Chirp 

Gini  Index 

0.0182 

0.2808 

0.4669 

0.6558 

0.6195 

0.6393 

Approved  for  Public  Release;  Distribution  Unlimited. 
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Table  5  Sparseness  Measure  of  Various  1-D  signals  (signal  length  =  1024). 


FFT 

DCT 

WAV 

ELET 

Haar 

Coiflet 

Db4 

Symiet 

Lq 

0.5918 

0.0176 

0.0449 

0.0146 

0.0186 

0.0156 

Lp  (p=l) 

10.6070 

0.3577 

0.4440 

0.2852 

0.3325 

0.3063 

Lp  (p=0.5) 

3506.9000 

98.7320 

124.6000 

8.0445 

11.0660 

8.7073 

Log 

0.8779 

0.0448 

0.0752 

0.0343 

0.0434 

0.0365 

Tanh 

0.9173 

0.0689 

0.II22 

0.0360 

0.0414 

0.0358 

Kurtosis 

0.3519 

0.2759 

0.3764 

0.I92I 

0.1582 

0.1592 

Heavisine 

Gini  Index 

0.8453 

0.8839 

0.8956 

0.9890 

0.9863 

0.9884 

0.6650 

0.0752 

0.0469 

0.0361 

0.0557 

0.0508 

Lp  (P=l) 

11.2080 

0.3167 

0.2098 

0.1926 

0.1958 

0.1936 

Lp  (p=0.5) 

7381.3 

189.07 

49.617 

33.173 

33.605 

32.931 

Log 

1.4435 

0.0761 

0.0583 

0.0537 

0.0552 

0.0542 

Tanh 

0.9815 

0.2041 

0.II30 

0.1064 

0.II09 

0.1089 

Kurtosis 

0.0327 

0.0390 

0.0440 

0.0512 

0.0472 

0.0500 

Bumps 

Gini  Index 

0.6603 

0.7041 

0.8971 

0.9158 

0.9160 

0.9160 

Lo 

0.7197 

0.1065 

0.0361 

0.0645 

0.0615 

0.0615 

Lp  (P=l) 

22.7540 

0.6584 

0.3909 

0.4747 

0.4812 

0.4694 

Lp  (p=0.5) 

I356I. 

353.67 

20.58 

67.345 

59.926 

65.456 

Log 

1. 9510 

0.1358 

0.0798 

0.1045 

0.1030 

0.1032 

Tanh 

0.9994 

0.3007 

0.0662 

0.1386 

0.1347 

0.1394 

Kurtosis 

0.1880 

0.2055 

0.1893 

0.1982 

0.1824 

0.2077 

Blocks 

Gini  Index 

0.6935 

0.7385 

0.9687 

0.9345 

0.9403 

0.9346 

Lo 

0.3604 

0.0146 

0.0166 

0.0098 

Lp  (P=l) 

2.9376 

0.0819 

0.0767 

0.0445 

0.0498 

0.0438 

Lp  (p=0.5) 

1405.8 

34.259 

30.87 

2.9673 

4.8014 

3.0399 

Log 

0.4948 

0.0155 

0.0147 

0.0I2I 

0.0125 

0.0II7 

Tanh 

0.6310 

0.0449 

0.0441 

0.0247 

0.0351 

0.0246 

Kurtosis 

0.0537 

0.0841 

0.1063 

0.1080 

0.I55I 

0.1305 

Doppler 

Gini  Index 

0.7862 

0.8140 

0.8324 

0.9717 

0.9593 

0.9710 

Lo 

0.4365 

0.0088 

0.0039 

0.0078 

Lp  (P=l) 

2.3767 

0.0698 

0.0359 

0.0337 

■iIiBIm 

0.0344 

Lp  (p=0.5) 

1366.5 

35.802 

5.6859 

0.9825 

0.7052 

0.9895 

Log 

0.4163 

0.0II9 

0.0087 

0.0091 

0.0086 

0.0095 

Tanh 

0.9316 

0.0262 

0.0153 

0.0175 

0.0144 

0.0179 

Kurtosis 

0.1667 

0.1507 

0.2470 

0.2336 

0.2769 

0.2598 

Ramp 

Gini  Index 

0.6878 

0.7420 

0.9483 

0.9862 

0.9889 

0.9859 

Lo 

0.7432 

0.3330 

0.2178 

0.1748 

0.1982 

0.I7I9 

Lp  (P=l) 

22.3430 

0.6276 

0.5285 

0.4239 

0.4396 

0.4285 

Lp  (p=0.5) 

22861.0 

586.11 

426.32 

235.48 

268.07 

241.42 

Log 

2.6977 

0.1602 

0.1373 

0.II96 

0.I2I0 

0.1208 

Tanh 

1. 0000 

0.7224 

0.5I7I 

0.3687 

0.3939 

0.3730 

Linear 

Kurtosis 

0.0010 

0.0015 

0.0035 

0.0056 

0.0074 

0.0055 

Chirp 

Gini  Index 

0.0135 

0.2787 

0.4732 

0.6576 

0.6243 

0.6504 

Approved  for  Public  Release;  Distribution  Unlimited. 
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0.22923 


0.24258 


0.0342 

0.1634 

61.465 

0.0426 

0.0790 

0.0431 

0.9259 


0.2442 


0.3862 


Table  6  Sparseness  Measures  of  1-D  Signals  (Signal  length  =2048,  sample=256). 


ymie 


Lo 

Lp  (p=l) 

L„(p=0.5) 

Log 

tanh  (a=2,  b=3) 
Kurtosis 

Gini  Index _ 

MSE-BP _ 

Heavisine  mse-omp 


Lfl 

L„(p=l) 

L„(p=0.5) 

Log 

tanh  (a=2,  b=3) 
Kurtosis 

Gini  Index _ 

MSE-BP _ 

Bumps  MSE-OMP 


Lfl 

Ln  (P=l) 
L„(p=0.5) 

Log 

tanh  (a=2,  b=3) 
Kurtosis 
Gini  Index 


MSE-BP 


MSE-OMP 


Lfl 

Ln  (P=l) 

L„(p=0.5) 

Log 

tanh  (a=2,  b=3) 
Kurtosis 

Gini  Index _ 

MSE-BP 

Doppler  MSE-OMP 


Lo 

Ln  (P=l) 
Ln(p=0.5) 

Log 

tanh  (a=2,  b=3) 
Kurtosis 
Gini  Index 


MSE-BP 


MSE-OMP 


Lfl 

Ln(P=l) 

Ln(p=0.5) 

Log 

tanh  (a=2,  b=3) 
Kurtosis 

Gini  Index _ 

MSE-BP _ 

Linear  Chirp  mse-OMP 


FFT 

DCT 

0.5889 

11.1700 

6903.00 

133.98 

0.8748 

0.0296 

0.9170 

0.0477 

0.3520 

0.2757 

0.8540 

0.8913 

9.5246 

0.04213 

0.013563 

0.083194 

0.6226 

0.0479 

11.9480 

0.2382 

12410.0 

224.95 

1.2137 

0.0570 

0.9479 

0.1321 

0.0325 

0.0388 

0.7648 

0.7927 

0.5198 

0.3227 

0.1329 

0.5740 

0.7256 

0.0728 

25.1980 

0.5197 

27473.0 

521.46 

1.9378 

0.0979 

0.9995 

0.2155 

0.1875 

0.2050 

0.7263 

0.7576 

6.0433 

1.0951 

0.22351 

1.0854 

0.2217 

0.0117 

2.9407 

0.0582 

1573.3 

28.841 

0.3922 

0.0118 

0.3679 

0.0316 

0.0537 

0.0841 

0.8914 

0.9027 

0.085698 

0.017692 

0.0040537 

0.030266 

0.0374 

0.1887 

0.9833 


0.21728 


1.02E-29 


0.4390 

2.5971 

2768.0 

0.4171 

0.9316 

0.1667 

0.7142 


0.099817 


0.0014906 


0.0005 

0.0098 


0.4914 


0.3728 


0.0183 

0.1502 

0.7638 


0.0052462 


0.010011 


.3330 
.6302 
1177.90 
0.1612 
0.7250 
0.0007 
0.2771 


0.6083 


1.1850 


.1065 

.8763 


0.014163 


0.021305 


0.0049 

0.0259 

6.0158 

0.0060 

0.0102 

0.2504 

0.9641 


0.00017928 


1.41E-04 


0.2227 

0.5267 

845.56 

0.1370 

0.5140 


0.0088 

0.0107 

0.0107 

0.2005 

0.2359 

0.2167 

6.0631 

8.2454 

6.5401 

0.0209 

0.0265 

0.0223 

0.0218 

0.0242 

0.0215 

0.1936 

0.1583 

0.1626 

0.9942 

0.9930 

0.9940 

0.005685 

0.0071205 

0.0050828 

8.73E-07 

2.61E-06 

1.78E-07 

0.0239 

0.0322 

0.0337 

0.1363 

0.1430 

0.1364 

26.204 

28.434 

26.146 

0.0367 

0.0391 

0.0366 

0.0633 

0.0684 

0.0638 

0.0537 

0.0468 

0.0524 

0.9554 

0.9538 

0.9555 

0.1867 

0.1686 

0.1842 

0.3495 

0.4997 

0.2243 

0.0391 

0.0386 

0.0391 

0.3401 

0.3485 

0.3416 

56.216 

50.545 

56.939 

0.0674 

0.0680 

0.0690 

0.0846 

0.0823 

0.0857 

0.1966 

0.1819 

0.2064 

0.9638 

0.9664 

0.9622 

0.43873 

0.63497 

0.46815 

0.38539 

0.40262 

0.48042 

0.0073 

0.0098 

0.0073 

0.0328 

0.0354 

0.0320 

2.4776 

3.8096 

2.4683 

0.0089 

0.0093 

0.0086 

0.0171 

0.0210 

0.0165 

0.1025 

0.1551 

0.1225 

0.9845 

0.9787 

0.9844 

0.00070946 

0.0027596 

0.00094719 

1.92E-05 

6.80E-04 

8.70E-05 

.0039 
.0235 
.6948 
.0061 
.0096 
0.2336 
0.9932 


1.21E-05 


1.87E-11 


0.5527 


1.1756 


0.5452 


1.1411 


.0034 
.0221 
.5388 
.0058 
0.0087 
0.2769 
0.9943 


2.29E-07 


2.90E-14 


0.1914 

0.4420 

532.37 

0.1227 

0.3925 

0.0030 

0.6276 


0.5486 


1.0732 


I 


.0039 
.0241 
.7179 
.0064 
0.0102 
0.2594 
0.9930 


2.00E-05 


1.67E-11 


.1670 
.4290 
77.22 
.1212 
0.3731 
0.0028 
0.6522 


0.5575 


1.2488 
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Table  7  Sparseness  IVI 

[easures  of 

l-D  Signals  (Signal  length  =4096,  sample=512). 

FFT 

DCT 

Haar 

by  m  let 

^0 

0.5933 

■HR 

0.0063 

Lp  (p=l) 

11.7330 

0.2292 

WBEm 

0.1536 

L„  (p=0.5) 

13650 

184.39 

146.36 

WBBm 

6.1765 

4.9255 

Log 

0.8732 

mSSM 

0.0321 

0.0136 

tanh  (a=2,  b=3) 

0.9175 

0.0450 

0.0127 

0.0128 

Kurtosis 

0.3519 

0.3760 

0.1941 

^R^^R 

0.1643 

Gini Index 

0.8614 

0.9476 

0.9970 

0.9969 

MSE-BP 

9.5268 

0.019721 

0.053432 

1.57E-05 

8.52E-06 

1.62E-05 

Heavisine 

MSE-OMP 

0.0075 

0.0187 

0.0321 

2.25E-06 

1.42E-06 

1.88E-06 

Lq 

0.4763 

0.0332 

0.0208 

0.0149 

0.0188 

0.0193 

(P=l) 

12.1710 

0.1715 

0.1223 

0.0998 

0.1016 

0.0988 

L„  (p=0.5) 

18003 

230.590 

71.547 

21.516 

21.879 

20.813 

Log 

0.9151 

0.0409 

0.0295 

0.0257 

0.0264 

0.0254 

tanh  (a=2,  b=3) 

0.8296 

0.0838 

0.0504 

0.0384 

0.0391 

0.0378 

Kurtosis 

0.0326 

0.0389 

0.0429 

0.0525 

0.0468 

0.0519 

Gini Index 

0.8545 

0.8708 

0.9477 

0.9757 

0.9761 

0.9761 

MSE-BP 

0.5196 

0.2010 

0.0593 

0.0804 

0.0470 

0.0794 

Bumps 

MSE-OMP 

0.0606 

0.0187 

0.0346 

0.0141 

0.0122 

0.0244 

0.7173 

0.0535 

0.0247 

(P=l) 

27.6350 

0.4018 

0.2514 

L„  (p=0.5) 

55271.0000 

736.4300 

^R^hR 

41.1600 

Log 

1.9268 

0.0698 

RRIgH 

0.0444 

RRQRR 

tanh  (a=2,  b=3) 

0.9954 

0.1519 

0.0506 

Kurtosis 

0.1872 

0.2047 

0.1817 

Gini Index 

0.7537 

0.7841 

0.9815 

MSE-BP 

6.055 

0.40919 

7.60E-09 

0.15354 

0.17859 

0.1532 

Blocks 

MSE-OMP 

0.1446 

0.4088 

1.41E-29 

0.0824 

0.0462 

0.0655 

0.1274 

0.0076 

0.0039 

0.0051 

0.0044 

(P=l) 

2.9415 

0.0419 

0.0231 

0.0250 

0.0228 

L„  (p=0.5) 

1688.6000 

22.6940 

42.3990 

1.7667 

2.8356 

1.7841 

Log 

0.2892 

■sn 

0.0061 

0.0066 

0.0059 

tanh  (a=2,  b=3) 

0.2053 

WBBm 

0.0098 

0.0120 

0.0100 

Kurtosis 

0.0537 

WBBM 

0.1551 

0.1186 

Gini Index 

0.9453 

0.9923 

0.9893 

0.9922 

MSE-BP 

0.085693 

0.006159 

0.0071591 

6.46E-05 

0.00027183 

9.24E-05 

Doppler 

MSE-OMP 

0.0018 

0.0065 

0.0053 

5.06E-06 

3.33E-05 

5.16E-06 

Lq 

0.4387 

0.0042 

0.0027 

■SEH 

0.0017 

0.0024 

(P=l) 

2.8176 

0.0411 

0.0185 

0.0158 

0.0172 

L„  (p=0.5) 

5585.6000 

72.4820 

6.2066 

0.4087 

0.5595 

Log 

0.4175 

0.0061 

0.0040 

■■R 

0.0038 

0.0043 

tanh  (a=2,  b=3) 

0.9316 

0.0130 

0.0061 

0.0060 

0.0050 

0.0065 

Kurtosis 

0.1667 

0.1500 

0.2522 

0.2337 

0.2769 

0.2593 

Gini  Index 

0.7365 

0.7819 

0.9750 

0.9965 

0.9970 

0.9963 

MSE-BP 

0.10012 

0.0028821 

6.78E-05 

8.58E-09 

2.25E-09 

3.39E-09 

Ramp 

MSE-OMP 

0.0009 

0.0032 

2.54E-05 

1.28E-11 

4.82E-14 

9.97E-12 

0.7466 

0.3276 

0.2188 

0.1694 

0.1704 

(P=l) 

44.9710 

0.6321 

0.5263 

0.4229 

0.4269 

L„  (p=0.5) 

184120.0000 

2363.5000 

1683.5000 

915.4600 

940.0700 

Log 

3.3053 

0.1619 

0.1371 

0.1200 

RSIhRI 

0.1207 

tanh  (a=2,  b=3) 

1.0000 

0.7268 

0.5102 

0.3662 

0.3700 

Kurtosis 

0.0002 

0.0004 

0.0010 

0.0015 

0.0016 

Gini  Index 

0.0071 

0.2761 

0.4801 

0.6624 

0.6555 

MSE-BP 

0.4943 

0.5862 

0.5781 

0.6071 

0.5748 

0.5872 

Linear  Chirp 

MSE-OMP 

0.4044 

0.9545 

1.0680 

0.9360 

0.9230 

1.0400 
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The  reconstructed  block  and  Heavisine  signal  from  an  original  signal  length  of 
4096,  using  400  CS  samples  are  shown  in  Fig:ure(s)  8-9. 


Original  Signal 


FFT 


Coiflet  Wavelet 


Figure  8  Block  signal  reconstruction,  signal  length  4096,  number  of  samples  400. 


Original  Signal 


0  1000  2000  3000  4000 


Haar  Wavelet 


Coiflet  Wavelet 


0  1000  2000  3000  4000 

Symiet  Wavelet 


Figure  9  Heavisine  signal  reconstruction,  signal  length  4096,  number  of  samples  400. 

2.2.3  Results  and  Discussions  of  1-D  Compressive  Sampling  Signal  Reconstruction 
with  Noise 


For  the  noisy  situation,  two  different  scenarios  were  considered  and  they  are 
illustrated  in  Figure(s)  10  and  11.  In  the  first  scenario  shown  in  Figure  10,  an  additive 
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white  Gaussian  noisy  channel  was  considered.  In  this  case,  it  was  assumed  that  the 
original  signal  does  not  contain  noise  and  noise  was  added  to  the  decimated  data  during 
transmission.  Hence  the  result  of  the  sparseness  measure  is  not  affected.  However,  we 
tried  to  investigate  which  basis  result  in  the  best  signal  reconstruction,  even  after  the 


transmitted  decimated  samples  have  been  corrupted  by  noise. 


Figure  10  CS  Signal  reconstruction  in  a  Noisy  Channel. 


In  Figure  II,  a  situation  whereby  the  original  signal  is  corrupted  by  Gaussian  noise 
before  being  transformed  by  the  basis  was  considered.  In  this  case,  the  sparseness 
measures  were  re-evaluated  and  the  results  obtained  for  the  10  dB  noise  level  in  the 
original  signal  for  the  Block  signal,  are  shown  in  Figure(s)  12-  16.  All  the  sparseness 
measures  predicted  that  the  best  basis  for  the  block  signal  reconstruction  is  Haar  wavelet, 
only  the  Kurtosis  measure  predicted  otherwise,  the  plot  from  the  kurtosis  was  not 
included. 
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Ljj  norm  of  Block  Signal  (Noisy)  vs  Signal  Length 


Figure  12  Blocks  Sparsity  -  lo  Pseudo-Norm  (noisy). 


norm  of  Block  Signal  (Noisy)  vs  Signal  Length 


Figure  13  Blocks  Sparsity  -  h  Norm  (noisy). 
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Log  Measure 


Log  Measure  of  Block  Signal  (Noisy)  vs  Signal  Length 


Figure  14  Blocks  Sparsity  -  Log  Measure  (noisy). 


Tanh  of  Block  Signal  (Noisy)  vs  Signal  Length 


Figure  15  Blocks  Sparsity  -  Tanh  Measure  (noisy). 
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Gini  Index  of  Block  Signal  (Noisy)  vs  Signal  Length 


Figure  16  Blocks  Sparsity  -  Gini  Index  Measure. 


Original  Signal  with  noise 

FFT 

y 

ftnwifi  I  \T\ 
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4 

2 

0 

-2 

DOT 

Coiflet  Wavelet 

6 

4 
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0 

-2 

Db4  Wavelet 

Symlet  Wavelet 

6 

4 

0 

-2 

Haar  Wavelet 

original  signal  without  noise 

6 

4 

2 

Ah 

; 

0 

-2 

L 

Figure  17  Noisy  Block  signal  reconstruction,  signal  length  4096,  samples  size  400. 

For  each  one  of  the  sparse  signals  considered  earlier,  the  effect  of  a  noisy 
communication  channel  on  the  reconstructed  signals  using  compressive  sampling  based 
signal  reconstruction  was  also  investigated.  At  different  noise  level,  using  the  orthogonal 
matching  pursuit  algorithm,  the  results  obtained  were  shown  in  Figure(s)  18  to  22.  For 
the  Heavisine  signal  type  at  different  SNR,  the  bases  with  the  lowest  MSB  are  symlet. 

Approved  for  Public  Release;  Distribution  Unlimited. 


23 


coiflet  and  Db4  which  correlates  with  the  earlier  prediction  of  the  sparseness  measures. 
Figure  19  is  the  plot  obtained  for  the  bump  signal  while  the  plot  for  the  block  signal  is 
shown  in  Figure  20.  For  the  block  signal,  under  different  noise  condition,  the  haar 
wavelet  was  clearly  the  best  basis  for  signal  reconstruction  as  predicted  by  the  sparseness 
measures  while  for  the  bump  signal,  it  can  be  seen  that  the  DCT  is  obviously  not  a  good 
sparsifying  basis  to  be  used.  Similar  results  were  obtained  for  the  Doppler  and  ramp 
signals  and  the  plots  are  shown  in  Figure(s)  21  and  22  respectively. 


OMP  Heavisine  Signal  mse  vs  noise  level 


Figure  18  Heavisine  Signal-mean  square  error  vs  SNR  using  OMP. 
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Mean  Square  Error  Mean  Square  Error 


OMP  Bumps  Signal  mse  vs  noise  level 


Figure  19  Bumps  Signal-mean  square  error  vs  SNR  using  OMP. 


OMP  Blocks  Signal  mse  vs  noise  level 


Figure  20  Blocks  Signal  mean  square  error  vs  SNR  using  OMP. 
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Mean  Square  Error  Mean  Square  Error 


OMP  Doppler  Signal  mse  vs  noise  level 


Figure  21  Doppler  Signal-mean  square  error  vs  SNR  using  OMP. 


OMP  Ramp  Signal  mse  vs  noise  level 


Figure  22  Ramp  Signal  mean  square  error  vs  SNR  using  OMP. 
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Observations  from  the  various  results  and  data  obtained  from  the  simulation  show  that 
although  the  CS  measurements  may  be  performed  with  no  a  priori  knowledge  about  the 
signal  (other  than  the  fact  that  the  signal  is  sparse  or  compressible  in  some  basis),  the 
reconstruction  error  can  depend  strongly  on  the  choice  of  basis  functions.  Hence  a  means 
for  determining  the  best  sparsifying  basis  function  is  necessary.  It  has  been  shown  that  it 
is  possible  to  predict  which  sparsifying  basis  will  result  into  the  best  CS  signal 
reconstruction.  With  the  result  obtained  from  all  the  signals,  it  was  discovered  that  all  of 
the  sparsity  measures  performed  reasonably  well  except  for  kurtosis,  with  lo  norm  and 
Gini  index  giving  a  much  better  performance.  The  Gini  index  has  some  advantages  over 
the  other  measures.  In  addition  to  providing  a  good  indication  of  signal  sparsity,  it  does 
not  require  normalization.  Its  value  always  lies  between  0  and  1 ,  regardless  of  the  signal 
length,  and  is  thus  a  better  pointer  for  decision  making  when  different  signal  lengths  are 
being  considered.  The  lo  norm  measure,  apart  from  requiring  normalization,  also  requires 
rounding  off  the  values,  since  most  of  the  data  after  transformation  are  not  truly  zero. 

2,3  Nyquist  folding  analog-to-information  receiver  in  fading  environment 

In  Cognitive  Radio  applications  as  well  as  battlefield  communication  systems,  dynamic 
spectrum  access  or  detection  is  of  great  importance.  In  the  United  States,  the  FCC  recently 
made  a  decision  to  allow  white  spaces  to  be  used  by  unlicensed  devices.  This  mling  will 
make  available,  the  frequency  band  formally  being  used  by  TV  stations  for  analog 
transmission.  The  ability  to  elficiently  and  effectively,  sense  or  detect  the  presence  of 
transmitted  signal  at  specific  frequencies  and  bandwidth  in  the  wide  banded  spectrum  will 
help  a  great  deal  in  each  of  the  application  areas  mentioned. 

Wideband  spectrum  monitoring  is  a  major  challenge  in  Cognitive  Radio  [22].  Current 
techniques  used  in  identifying  or  detecting  white  spaces  in  the  frequency  spectrum  include; 
transmitter  detection,  cooperative  detection,  and  interference-based  detection  [23]. 
Transmitter  detection  approach  focuses  on  detecting  weak  signals  from  primary  transmitter 
by  carrying  out  observation  at  the  end  user.  There  exist  three  main  techniques  of  achieving 
this,  namely;  matched  filter  detection,  energy  detection  [24],  and  cyclostationarity  feature 
detection  [25]. 

Opportunistic  transmission  and  detection  requires  the  ability  to  detect  holes  in  the 
spectrum,  effectively  and  efficiently  makes  use  of  the  available  spectrum  by  sending  as 
many  information  as  possible  within  the  shortest  period  of  time.  In  the  past  few  decades, 
there  has  been  tremendous  research  efforts  in  the  area  of  digital  signal  processing,  however, 
the  performance  of  analog-to-digital  converter  (ADC)  has  been  much  slower,  averaging 
roughly  1.5  bits  /  8  years  [26],  [27]. 

Compressive  Sampling  (CS)  can  be  used  to  collect  data  at  rates  lower  than  what  Nyquist 
specifies,  hence  can  help  in  overcoming  the  limitations  of  ADC  in  wideband  spectrum 
monitoring  for  CR  applications  [28].  The  underlying  principle  of  CS  is  that  signals  that  are 
known  to  be  sparse  in  some  known  basis,  can  be  sampled  based  on  the  information  rate  as 
opposed  to  the  Nyquist  criteria  of  twice  the  total  bandwidth.  In  this  case,  a  sample  rate 
much  lower  than  that  implemented  by  conventional  digital  receiver  techniques  is  possible  as 
long  as  the  spectrum  is  sparsely  occupied.  This  information  sampling  and  recovery  idea  is 
also  referred  to  as  analog-to-information  [27]  (A2I). 
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The  Nyquist  Folding  Receiver  (NYFR)  has  recently  been  proposed  as  an  A2I  approach 
for  sparse  wideband  RF  applications  [29].  While  many  CS  techniques  involve  random 
projections  followed  by  a  computationally  expensive  reconstruction  [28],  [l]-[2],  the 
NYFR  performs  analog  compression  via  a  non-uniform  sampling  process  that  induces  a 
chirp-like  modulation  on  each  received  signal.  This  induced  modulation  can  be  measured 
using  time-frequency  analysis  techniques  to  determine  the  original  RF  band  of  origin 
without  full  signal  reconstruction.  In  this  paper,  we  investigate  the  viability  of  achieving 
signal  detection  and  identification  of  original  RF  band  of  signals  present  in  the  NYFR 
data,  when  the  original  signal  passes  through  the  presence  of  multiplicative  (Rice)  and 
additive  (white  Gaussian)  noise. 

2.3.1  Nyquist  Folding  Analog-to-Information  Receiver 

Physical  realizations  of  CS  receivers  to  perform  A2I  conversion  include  random 
filters  [30],  random  pre-modulation  followed  by  integration  [28],  [31],  structured 
random  sampling  [32],  and  the  NYFR  [29]. 

The  NYFR  performs  pulsed  sampling  directly  at  RF  using  a  frequency  modulated 
carrier  sample  clock  and  allows  multiple  Nyquist  zones  (defined  by  half  of  the 
average  RF  sample  rate)  to  fold  into  an  analog  interpolation  filter  as  shown  in  Figure 
2.23.  Figure  2.23  illustrates  noise  being  added  to  the  original  signal  during 
transmission.  The  NYFR  does  not  use  an  anti-aliasing  filter  prior  to  sampling;  instead 
it  makes  use  of  a  wideband  pre-select  filter  H  (co).  The  interpolation  filter  makes  it 
possible  to  decouple  the  ADC  clock  from  the  sampling  clock.  While  the  ADC  clock 
is  uniform,  the  sampling  clock  is  not,  it  is  modulated. 


Figure  23  Nyquist  Folding  A-to-I  Receiver  Architecture 

The  RF  sample  times  correspond  to  the  zero  crossing  time  for  the  rising  voltage 
of  the  modulated  sample  clock.  As  shown  in  [29],  the  received  signal  has  an  induced 
modulation  of  the  same  form  as  the  RF  sample  clock  modulation  with  a  modulation 
scale  factor  that  depends  on  the  Nyquist  zone  of  the  signal.  In  the  special  case  of  a 
narrowband  modulated  (phase  or  frequency)  signal  (p  (t)  with  carrier  frequency  cosi 
such  that 

(pit)  =  (OJ  +  0{t)  (28) 

Supposing  the  input  signal  is  a  narrowband  signal  with  center  frequency  coc  and 
information  modulation  v|/  (t)  defined  by 

X  '(t)  =  cosicDj  +  y/it))  (29) 

Then  the  normalized  pre-envelope  at  the  interpolation  filter  output  will  be  [29] : 
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y,(t) 


j(ajt+ii/(t)-Me(t)) 

^  ? 


M  ^sgn{Ns-Np)kfj 

kn  —  \ \Np  -  //j  I  /  2~|  (|~*~|  denotes  eeil) 


(30) 


where  cOd  is  the  carrier  (folding  frequency)  after  translation  by  bandpass 
sampling,  M  is  the  modulation  scale  factor,  Icr  is  the  harmonic  selected  by  the 
interpolation  filter,  and  Np  and  NS  are  the  signed  Nyquist  zones  of  the  interpolation 
filter  and  RF  signal  respectively  relative  to  the  RF  sample  rate  fsi.  The  signed 
Nyquist  zone  for  a  signal  or  interpolation  filter  centered  at  fc  is  defined  by 


N  =  ±[2/p//^J  ([*J  denotes  floor)  (31) 

where  the  sign  is  negative  for  odd  (spectrally  reversed)  Nyquist  zones  and  even 
otherwise. 

Therefore,  though  multiple  signals  from  different  Nyquist  zones  may  alias  into  the 
same  band,  the  information  from  the  different  signals,  including  the  original  RF,  can 
still  be  recovered  based  on  the  fact  that  the  added  modulation  is  different  for  each 
Nyquist  zone  under  the  condition  that  the  signal  environment  is  relatively  sparse  so 
that  signals  are  separable.  Because  the  NYFR  does  not  rely  on  random  projections, 
the  NYFR  output  may  be  processed  using  time  frequency  methods  rather  than 
complex  1 1  minimization. 

In  order  to  analyze  the  data  obtained  from  the  NYFR  system,  a  time -frequency 
analysis  of  the  data  has  to  be  carried  out.  The  spectrogram  which  is  based  on  Fast 
Fourier  transform  (FFT)  has  been  widely  used  as  one  of  the  methods  for  time- 
varying  spectral  analysis.  It  has  been  applied  for  signal  analysis  in  many  applications 
such  as  radar,  sonar,  speech,  geophysics  and  biological  signals  [33].  The  wavelet 
analysis  which  is  well  suited  for  analysis  of  transient  signals  can  be  used  as  well. 
Wavelet  uses  time-frequency  atoms  of  different  sizes  that  depend  on  the  frequency  to 
break  down  the  signal  [34].  Both  FFT  spectrogram  and  AWT  scalogram  has  been 
applied  to  the  analysis  of  the  NYFR  data  [29],  [35]. 


2.3.2  Ricean  Noise  and  Gaussian  Noise 

The  earlier  work  [29]  did  not  consider  the  behavior  of  the  NYFR  under  noisy 
condition.  We  investigated  the  effect  of  Rice  noise  and  Gaussian  noise  on  the 
performance  of  the  NYFR.  Gaussian  noise  is  the  most  widely  used  noise  model  for 
communication  system  analysis.  The  Rice  noise  has  amplitude  which  is  Rice 
distributed.  Suppose  yi  and  yi  represents  two  Gaussian  variables,  with  common 
variance  and  mean  values  mj  and  m2,  the  non-centrality  parameter  v 
equals  mf  The  pdf  of  the  Rice  noise  is  modeled  as 


f{^)- 


—  1 


^^3 


2(t" 


x>0 


(32) 


0 


X  <  0 


Where  the  function  Io(u)  is  the  modified  zeroth  order  Bessel  function  of  the  first 
kind  defined  as 
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The  mean  value  of  the  Rieean  noise  is  (2+m)  and  the  varianee  is  4(l+2m).  The 
original  signal  was  passed  through  a  eombination  of  both  multiplieative  Rieean  noise 
and  additive  white  Gaussian  noise. 

2.3.3  Automatic  Detection  Algorithm 

FFT  spectrogram  plots  and  AWT  scalogram  plots  can  be  used  to  determine  both  the 
folding  frequency  and  the  induced  modulation  bandwidth.  This  often  requires  visual 
inspection  of  the  plotted  images.  For  the  purpose  of  automation  and  also  to  remove 
human  error  in  interpretation  of  the  plots,  we  worked  on  developing  an  algorithm  that  is 
capable  of  automatically  predicting  the  folding  frequency  of  the  signals  after  passing 
through  the  NYFR  system  as  well  as  the  induced  modulation  bandwidth.  The  algorithm 
worked  on  either  the  spectrogram  or  the  scalogram  data.  The  pseudocode  for  the 
algorithm  is  listed  below. 


Step  1 ;  Normalize  the  data  from  the  FFT  spectrogram  or  AWT  scalogram.  Let  us  assume 
that  the  data  is  A,  the  normalized  data  is  define  by 


(34) 


where  p>l.  The  effect  of  p  will  be  investigated  in  detail  in  later  publication. 

Step  2;  Compare  the  data  with  a  threshold  y*,  which  was  determined  heuristically.  This 
will  be  discussed  further  in  the  experimental  result  section. 

Step  3:  Set  values  below  the  threshold  to  zero  and  values  >  the  threshold  to  1.  Values  at 
the  detected  signal  points  should  be  greater  than  the  threshold  value. 

Step  4:  Compute  the  difference  matrix  of  the  new  data  and  find  the  location  of  the  non¬ 
zero  elements.  This  can  be  used  in  tracking  the  point  at  which  the  signal  begins  and 
ends. 


Step  5:  Determine  the  maximum  and  minimum  point  for  each  detected  signal. 

Step  6:  the  frequency  at  which  the  various  signals  folds  to,  the  mid-point  of  the  maximum 
and  minimum  point. 

Step  7;  Induced  modulation  bandwidth  is  the  difference  between  the  maximum  and 
minimum  points  for  each  signals. 

A  flow  chart  illustration  of  the  algorithm  is  shown  in  Figure  24. 
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Start 


Figure  24  Flow  Chart  for  the  detection  algorithm 

2.3.4  Results  and  Discussions  of  Nyquist  Folding  A-2-I  Receiver 
In  this  section,  we  examine  some  experimental  results  using  simulated  NYFR  data.  In 
this  example,  the  RF  sample  clock  ramps  up  and  down  between  1995  Msps  and  2005 
Msps  over  a  2  psec  period,  resulting  in  a  sample  rate  of  2  Gsps  +/-  10  MHz.  There  are 
two  signals  in  the  spectrum  originally,  a  10.3GHz  signal  as  well  as  a  10.6GHz  signal.  Fig. 
3  shows  the  FFT  spectrogram  of  the  NYFR  data  in  a  noiseless  situation. 

The  signal  at  10.3  GHz  will  fold  to  300  MHz  with  an  induced  modulation  bandwidth  of 
+/-  25  MHz  (M  =  round  (10.3  /  2.0)  =  5).  The  second  signal  at  10.6GHz  will  fold  to  600 
MHz,  with  an  induced  modulation  bandwidth  of  +/-  25  MHz  (M  =  round  (10.6  /  2.0)  =  5). 

The  earlier  work  did  not  consider  the  behavior  of  the  NYFR  under  multiplicative  noise 
condition.  We  investigated  the  effect  of  Ricean  noise  and  Gaussian  noise  on  the 
performance  of  the  NYFR.  Figure  25  shows  the  plot  obtained  for  a  situation  based  on 
Figure  23  (a  combination  of  Ricean  multiplicative  noise  and  additive  Gaussian  noise) 
whereby  the  SNR  is  3  dB  (SNR  was  measured  by  comparing  x’(t)  to  x(t)  in  Figure  23. 
Plots  obtained  under  only  Gaussian  or  Rice  channel  conditions  are  similar  to  Figure  26. 
We  carried  out  simulations  under  different  signal-to-noise  ratio  conditions.  Table  8  lists 
the  result  obtained  from  our  analysis. 
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FFT  Spectrogram  (Noiseless ) 


900 
800 
700 

1  600 

2 

“  500 
® 

D 

2  400 

li. 

300 
200 
100 
0 

0  1  2  3  4  5  6  7  8  9  10 

Time,  ps 

Figure  25  Spectrogram  of  data  (Noiseless) 
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Figure  26  Spectrogram  of  data  (Gaussian  and  Rice  channel) 
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Table  8:  Percentage  Error  Analysis 


Actual 

Carrier 

Freq. 

(GHz) 

Folding  Frequency 

Induced  Modulation  Bandwidth 

Estimated 

Carrier 

Freq. 

(GHz) 

Carrier 
Freq.  Est. 
Error  (%) 

SNR 

Theory 

MHz 

Estimated 

Est. 

Error  % 

Theory 

MHz 

Estimated 

Est. 

Error  % 

10.3 

Noiseless 

300 

301 

0.3 

50 

46.9 

6.2 

10.301 

0.01 

9dB 

300 

302 

0.6 

50 

46.7 

6.6 

10.302 

0.02 

6dB 

300 

304 

1.2 

50 

52.4 

4.8 

10.304 

0.04 

3  dB 

300 

305 

1.6 

50 

53.4 

6.8 

10.305 

0.05 

10.6 

Noiseless 

600 

599 

0.2 

50 

46.9 

6.2 

10.599 

0.01 

9dB 

600 

597 

0.5 

50 

46.8 

6.4 

10.597 

0.03 

6dB 

600 

605 

0.8 

50 

52.6 

5.2 

10.605 

0.05 

3  dB 

600 

607 

1.2 

50 

52.8 

5.6 

10.607 

0.07 

We  have  been  able  to  show  that  the  NYFR  system  is  robust  under  Rice  and  Gaussian 
noisy  channel  condition.  The  accuracy  of  the  prediction  is  very  high  and  the  automatic 
detection  algorithm  eliminates  error  associated  with  visual  analyzes  of  the  result  and  can 
aid  in  automated  and  embedded  system.  We  also  intend  to  look  at  detection  capability 
under  other  types  of  channel  conditions  such  as  Nakagami  fading  channel. 

2,4  Dynamic  Spectrum  Sensing  Using  Compressive  Sensing  Algorithms 

2.4.1  Problem  Statement 

Assuming  that  a  cognitive  radio  receiver  receives  a  signal  occupying  N  consecutive 
spectrum  band  with  frequency  boundaries  located  at  /o  </;<...  .  /„.  The  frequency 
response  of  the  received  signal  can  be  categorized  based  on  the  level  of  the  power 
spectral  density  of  the  received  signal  r(t)  and  is  illustrated  in  Figure  27.  If  the  PSD  level 
is  high,  medium  or  low;  the  spectrum  can  be  considered  as  black,  grey  or  white.  Whereby 
white  hole  and  sometimes  gray  spaces  can  be  opportunistically  used  by  the  CR  for 
transmission,  while  the  black  hole  are  not  to  be  used  in  order  to  avoid  interferences. 

In  a  wideband  frequency  spectrum  situation,  following  Nyquist  sampling  criteria,  very 
large  samples  will  be  required  to  effectively  recover  signal  r(t).  The  goal  of  each 
Cognitive  radio  is  to  effectively  sense,  classify  and  estimate  the  frequency  spectrum  of 
signal  r(t).  Our  focus  is  to  be  able  to  sense  and  classify  (identify  occupied  frequency 
bands)  the  signal  r(t). 


Approved  for  Public  Release;  Distribution  Unlimited. 


33 


"O 

C/D 

O 


Wideband  frequency  of  interest 


Figure  27  PSD  of  a  wideband  spectrum 


2.4.2  Results  and  Discussions  of  Dynamic  Spectrum  Sensing 

In  our  simulations  based  on  the  theory  of  compressive  sensing,  we  considered  a  wideband 
of  interest  with  the  signal  of  interest  being  the  frequency  domain  PSD  of  the  spectrum. 
We  employed  the  idea  of  CS  signal  reconstruction  in  order  to  detect  occupied/unoccupied 
spectral  bands.  We  revisited  the  problem  discussed  in  [36]  and  solve  the  problem  using 
OMP,  BP  and  LASSO  reconstruction  algorithm.  Our  simulations  revealed  that 
reconstruction  of  the  signal  is  possible  using  any  of  these  three  algorithms.  Other  CS 
algorithms  will  be  investigated  in  our  next  work 
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Figure  28  Noiseless  Signal  frequency  response  at  a  compression  ratio  of  70% 
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First  we  consider  a  noiseless  condition  and  investigated  reconstruction  at  different 
sampling  rates  using  the  earlier  mentioned  reconstruction  algorithms.  At  a  sampling  ratio 
of  70%,  the  result  obtained  from  the  three  reconstruction  algorithms  that  were  considered 
is  shown  in  Figure  27. 

We  also  investigated  the  effect  of  an  AWGN  channel  condition  on  the  detection  through 
CS  reconstruction  and  the  result  obtained  for  a  situation  whereby  the  signal  to  noise  ratio 
is  lOdB  is  shown  in  Figures  29. 
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Figure  29a  Noisy  Signal  (10  dB)  frequency  response  at  a  compression  ratio  of  70% 
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Figure  29b  Noisy  signal  (10  dB)  frequency  response  at  a  compression  ratio  of  80% 
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Figure  29c  Noisy  signal  (10  dB)  frequency  response  at  a  compression  ratio  of  20% 


Figures  30  shows  the  characterization  result  obtained  using  Basis  Pursuit  algorithm  under 
different  signal  to  noise  ratio  for  an  additive  white  Gaussian  noise  channel  condition.  It  is 
a  plot  of  the  root  mean  square  error  between  the  reconstructed  spectrum  and  the  original 
spectrum  at  different  noise  level  for  different  sampling  rates. 


Figure  30a  RMS  error  using  BP  algorithm  and  different  no.  of  samples  for  AWGN 

channel  at  1  iteration 
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Figure  30b  RMS  error  using  BP  algorithm  and  different  no.  of  samples  for  AWGN 

channel  at  200  iterations 


Figure  30c  RMS  error  using  BP  algorithm  and  different  lower  no.  of  samples  for 

AWGN  channel  at  5  iterations 
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Figure  30d  RMS  error  using  BP  algorithm  and  different  lower  no.  of  samples  for 

AWGN  channel  at  200  iterations 

Multipath  Fading  Channel 

A  principal  feature  of  a  cognitive  radio  (CR)  is  a  reliable  speetrum  sensing 
technique  whieh  enables  the  CR  to  detect  white  spaces  in  the  frequeney  band.  This  allows 
opportunistic  access  of  the  unlicensed  (secondary)  users  to  these  unoccupied  spaces 
without  causing  interference  to  licensed  (primary)  users.  In  many  scenarios  the  CR  may 
operate  in  a  multipath  fading  environment  where  speetrum  sensing  must  deal  with  the 
fading  effects  of  the  unknown  primary  PSD  signal.  It  important  to  study  the  effeets  of 
multipath  fading  on  the  performance  of  Compressive  Sensing  (CS)  based  speetrum 
sensing  algorithms. 

The  deteetion  performanee  of  a  CR  network  system  is  eritieally  affected  by  the 
multi-path  characteristic  of  the  environment.  In  a  high  speed  wideband  CR  network, 
multi -path  fading  is  one  of  the  most  serious  transmission  difficulties.  However,  there  are 
also  long  term  fluctuations  in  the  characteristics  of  the  wireless  channel  depending  on  the 
actual  position  of  the  Cognitive  Radio  and  the  mobile  station,  that  is,  in  urban  areas  the 
signal  is  affected  by  multi-path  propagation.  Theorieal  results  show  that  the  channel 
capacity  in  a  multi-path  fading  environment  is  always  lower  than  that  in  a  Gaussian  noise 
environment.  In  a  wideband  frequency  communication  channel  that  has  fading 
characteristies,  it  is  very  important  to  know  what  the  degradations  are  in  channel  capaeity 
due  to  multipath  fading,  and  also  to  what  degree  the  diversity  schemes  ean  bring  the 
channel  capacity  up  or  down.  The  previous  chapter  centered  on  the  effect  of  additive 
noise  to  CS  based  spectrum  detection.  This  chapter  is  centered  on  invetigating  the 
robustness  of  CS  algorithms  in  speetrum  sensing  in  the  presence  of  additive  and 
multplicative  noise. 

The  channel  amplitude  in  a  Rayleigh  fading  environment  is  always  lower  than 
that  of  a  Gaussian  noise  environment.  The  fading  ehannels  are  modeled  as  frequency 
non-selective  slow  Rayleigh  and  Riee  fading  channels  corrupted  by  AWGN.  Rayleigh 
and  Riee  fading  channels  are  useful  models  of  real-world  phenomena  in  wireless 
communications.  These  phenomena  include  multipath  scattering  effects,  time  dispersion. 


Approved  for  Public  Release;  Distribution  Unlimited. 


38 


and  Doppler  shifts  that  arise  from  relative  motion  between  the  transmitter  and  reeeiver. 
The  transmission  path  between  a  transmitter  and  reeeiver  ean  vary  from  a  simple  line-of- 
sight  (LOS)  path  to  a  non-line-of-sight  (NLOS)  path,  which  is  a  path  that  is  severely 
obstructed  by  buildings,  mountains,  and  foliage.  Furthermore,  when  a  CR  moves  in 
space,  the  speed  of  motion  impacts  how  rapidly  the  signal  level  fades. 

Two  widely  used  models  for  fading  channels  are  the  Rayleigh  fading  channel  and 
the  Rice  fading  channel,  where  the  Rayleigh  fading  channel  is  actually  a  special  case  of 
the  Rice  fading  channel.  The  Rayleigh  model  is  used  when  there  is  NLOS  between 
transmitter 

and  receiver,  and  all  of  the  received  signal  power  is  due  to  multipath.  The  Rice  model  is 
used  when  there  is  a  LOS  between  transmitter  and  receiver,  but  a  substantial  portion  of 
the  received  signal  power  is  also  due  to  multipath.  When  there  is  LOS  between 
transmitter  and  receiver  and  virtually  none  of  the  received  signal  power  is  due  to  multi- 
path,  the  non-fading  channel  model  is  used. 

Noise 


Fading 


Figure  31  CS  Based  Spectrum  Detection  in  an  AWGN  and  Fading  Channel 

Figure  31  illustrates  CS  based  spectrum  detection  in  an  AWGN  plus  fading 
channel.  The  process  involves  a  noiseless  PSD  passed  through  a  AWGN  +  fading 
channel.  The  noisy  PSD  signal  is  transformed  using  Haar  wavelet  sparsying  basis  and 
fewer  coefficients  are  taken.  Random  projection  of  the  signal  ensures  that  the  PSD  can  be 
perfectly  reconstructed  using  CS  reconstruction  algorithms. 

Rayleigh  Fading  Channel 

Rayleigh  fading  is  a  type  ofsmall-scale  fading.  If  the  multipath  reflective  paths 
are  large  in  number  and  are  all  NLOS,  the  envelop  of  the  received  signal  can  be 
statistcally  expressed  by  using  a  Rayleigh  fading  distribution.  The  Rayleigh  fading 
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distribution  has  a  probability  density  function  given  by 


p{r)  = 


f 


cjT. 


exp 


2cr^ 


>0 


(35) 


0,  otherwise 

Where  r  is  the  envelop  amplitude  of  the  received  signal,  and  a  is  the  time- 
average  power  of  the  multipath  signal.  The  probability  of  the  envelop  of  the  received 
signal  not  exceeding  a  specified  value  of  R,  is  given  by  the  expression 

P(R)  =Pr(r<R)  (36) 


=  i^V(T)dr 


(37) 


1- 


(38) 


The  Rayleigh  faded  component  is  also  known  as  random,  scatter,  or  diffuse 
component  [22], 

In  rural  areas  there  are  only  a  few  scatterers,  hence  the  channel  is  similar  to  an 
AWGN  channel. Therefore  the  performance  of  CS  algorithms  for  the  AWGN  channel  and 
the  Rayleigh  fading  channel  is  investigated.  A  Rayleigh-type  fading  with  arbitrary  fading 
rate  was  assumed  for  the  signal  to  get  an  understanding  of  the  system  performance  in  a 
fast  fading  environment.  Rayleigh  fading  occurs  when  a  transmitted  signal  passes  through 
a  medium  where  there  are  a  large  number  of  scatterers.  The  results  (Figures  32  and  33) 
are  generated  to  show  the  degradation  of  channel  capacity  in  a  Rayleigh  fading  and 
additive  noise  environment  and  the  performance  of  CS  algorithms  for  spectrum  hole 
detection. 
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Figure  32  Noisy  PSD  (Rayleigh)  at  a  Compression  Ratio  of  60%. 


Figure  33  Noisy  PSD  (Rayleigh)  at  a  Compression  Ratio  of  80%. 
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Figure  34  RMS  Error  Using  CS  Algorithms  of  Noisy  PSD  (10  dB)  at  70%  Sampling 

Measurement 


Figures  32  and  33  shows  the  performance  of  BP,  OMP,  and  LASSO  in  a  noisy 
AWGN  and  Rayleigh  fading  channel.  The  effect  of  fading  rate  on  a  CR  network  system 
performance  was  evaluated.  Figure  33  gives  precise  information  on  the  amplitude 
degradation  of  the  reconstructed  spectrum  by  BP,  OMP,  and  LASSO  reconstruction 
algorithm.  This  amplitude  degradation  is  not  necessarily  a  problem  in  spectrum  sensing 
since  the  frequency  occupancy  is  distinctively  known.  Although  significant  error  occurs 
in  the  magnitude  of  the  characterization  based  on  Rayleigh  fading  channel  plus  AWGN 
in  comparison  to  just  AWGN  channel,  this  is  not  necessarily  a  problem  since  the  goal  is 
to  determine  the  occupied  frequency  band  in  the  spectrum. 

Root  mean  square  error  is  employed  in  Figure  34  as  a  performance  metric  to 
investigate  the  algorithm  with  the  best  performance.  Despite  the  fact  that  OMP  performed 
better  at  70%  sampling  ratio  the  three  algorithms  had  the  same  result  at  10  dB  SNR. 
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Figure  35  RMS  Error  Using  BP  Algorithm  at  Different  Number  of  Samples  for 

Noisy  Channel 


Figure  36  RMS  Error  Using  LASSO  Algorithm  at  Different  Number  of  Samples  for 

Noisy  Channel 
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Figure  37  RMS  Error  Using  OMP  Algorithm  at  Different  Number  of  Samples  for 

Noisy  Channel 

Figures  35  -  37  show  how  different  sample  sizes  affect  the  performance  of  the 
BP,  OMP,  and  LASSO  algorithm.  From  Figure  37,  it  can  be  deduced  that  OMP  performs 
better  with  sampling  size  above  70%  as  compared  to  BP  and  LASSO.  It  is  also  noted  that 
as  the  SNR  goes  from  0  to  20  dB,  the  effect  of  sample  size  on  the  reconstruction  of  the 
PSD  becomes  negligible  as  the  RMSE  converges.  Least  Absolute  Shrinkage  Selective 
Operator  performs  better  when  less  than  70%  compression  ratio  is  used.  The  results 
obtained  show  that  in  noisy  and  Rayleigh  fading  channels  the  CS  algorithms  have 
different  performances.  Depending  on  the  environment  of  CR  network  the  primary  user 
spectrum  occupancy  can  be  known  and  false  alarm  can  be  prevented. 

Rice  Fading  Channel 

Wireless  communication  channels  experience  fading.  The  fading  is  caused  by 
reflections,  diffractions  and  scatterings  in  the  channel.  The  small-scale  fading  is  Rice 
fading  when  there  is  a  dominant  non  fading  signal  component  present,  such  as  a  LOS 
propagation  path. 

The  Rice  fading  distribution  has  a  probability  density  function  given  by 

p(r)  =  exp  ^0  (J) ,  r  >  0  and  d  >  0  ;r  <  0  (39) 

Where  r  is  the  envelope  and  a  is  average  of  the  Rice  fading  signal.  A  is  the  peak 
amplitude  of  the  dominant  signal,  and  /o(.)  is  the  modified  Bessel  function  of  zero  order, 
which  is  given  by  [23] 

^0(Z)  ~  ^T1=0  22«n!n! 

Note  that  the  Rice  fading  distribution  in  equation  (39)  presents  two  extreme  cases: 
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1.  If  the  dominant  signal  A  =  0,  p(r)  become  the  Rayleigh  fading  distribution  in 
equation  (35). 

2.  If  the  dominant  signal  A  is  large,  p(r)  becomes  a  Gaussian  distribution. 

The  Rice  fading  distribution  is  usually  expressed  in  terms  of  a  parameter  K  that  is 
defined  as  the  ratio  between  the  deterministic  signal  power  and  the  variance  of  the 
multipath.  The  parameter  is  given  by 


k  = 


2<J' 


Or  in  decibels 


'^  =  201og(^). 


(41) 


(42) 


Where  k  (or  K)  is  known  as  the  Rice  factor  that  completely  specifies  the  Rice  fading 
distribution.  Note  that  the  Rice  fading  distribution  becomes  the  Rayleigh  fading 
distribution  as>l  K  —co  dB. 

The  Multipath  Rice  Fading  Channel  in  MatLab  is  used  to  model  mobile  wireless 
communication  systems  when  the  transmitted  signal  can  travel  to  the  receiver  along  a 
dominant  LOS  or  direct  path.  The  effects  of  an  AWGN  and  Rice  fading  channels  on  CS 
based  spectrum  detection  was  highlighted  and  the  result  obtained  for  a  situation  whereby 
the  signal  to  noise  ratio  is  lOdB  and  Rice  fading  K-factor  is  10  seconds  is  shown  in 
Figures  38  to  39. 

Comparing  the  result  of  Figure  40  and  41  shows  that  OMP  performs  better  at  a 
higher  signal  sampling  measurement,  but  the  performance  deteriorates  greatly  at  a 
sampling  measurement  below  70%  .This  characteristic  performance  of  OMP  can  be 
justified  theoretically  by  the  fact  that  OMP  is  a  compressed  sensing  solver  that  seeks  to 
compute  the  minimum  norm  solution. 
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Figure  38  Noisy  PSD  (Rice  fading  +  Gaussian)  at  a  Compression  Ratio  of  70%. 
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PSD  PSD 


Original  PSD 
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Figure  39  Noisy  PSD  (Rayleigh  Fading  Channel  +  Gaussian)  Frequency  Responses 

at  a  Compression  Ratio  of  60%. 


Signal  to  Noise  Ratio(dB) 


Figure  40  RMS  Error  Using  CS  Algorithms  of  Noisy  (Rice  Fading+AWGN)  PSD  (10 

dB)  at  70%  Sampling  Measurement 
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Figure  41  RMS  Error  Using  CS  Algorithms  of  Noisy  (Rice  fading+AWGN)  PSD  (10 

dB)  at  60%  Sampling  Measurement 

A  test  of  the  performance  of  CS  reconstruction  algorithms  for  spectrum  sensing  in 
a  Rice  fading  environment  is  accomplished  by  passing  the  PSD  through  Rice  fading 
channel  in  AWGN  channel.  The  reconstructed  PSD  is  as  shown  in  Figures  2.40  and  2.41. 


Figure  42  RMS  Error  Using  BP  algorithm  at  Different  Number  of  Samples  for  Rice 

Fading  plus  AWGN  Noise 
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Figure  43  RMS  Error  Using  BP  Algorithm  at  Different  Sample  Number  for  Rice 

Fading  plus  AWGN  Noise 


Figure  44  RMS  Error  Using  OMP  Algorithm  at  Different  Number  of  Samples  for 

Rice  Fading  plus  AWGN  Noise 
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Figure  45  RMS  Error  Using  OMP  Algorithm  at  Different  Sample  Number  for  Rice 

Fading  plus  AWGN  Noise 


Figure  46  RMS  Error  Using  UASSO  Algorithm  at  Different  Number  of  Samples  for 

Rice  fading  plus  AWGN  Noise 
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Figure  47  RMS  Error  Using  LASSO  Algorithm  at  Different  Sample  Number  for 

Rice  fading  plus  AWGN  Noise 


The  characterization  result  of  the  RMSE  using  BP,  OMP  and  LASSO  algorithms 
evaluated  at  different  sampling  rates  is  shown  in  Figures  40  to  47.  Due  to  the  random 
nature  of  CS  sampling,  the  recovery  procedure  is  done  one  hundred  times  for  each 
sample  size  of  every  test.  The  reconstruction  is  carried  out  using  BP,  OMP  and  LASSO. 
Root  mean  square  error  of  the  reconstructed  PSD  in  comparison  to  the  original  PSD  were 
computed  and  used  as  performance  metric. 

The  results  obtained  validate  the  fact  that  the  CS  approach  to  spectrum  detection  is 
robust  whether  in  an  AWGN  channel  or  a  Rice  fading  channel.  Although  considerable 
errors  occurs  in  the  magnitude  of  the  characterization  based  on  Rice  fading  channel  plus 
AWGN  in  comparison  to  AWGN  channel,  this  is  not  necessarily  a  problem  since  the 
goal  is  to  determine  the  unoccupied  frequency  band  in  the  spectrum. 

Doppler  Spread 

Due  to  the  relative  motion  of  transmitter  and  receiver  in  the  CR  networks,  the 
Doppler  shift  is  of  great  practical  importance  to  take  account  of  The  effect  of  the 
Doppler  shift  on  the  CR  network  receiver  is  prevalent.  In  a  dense  CR  network  the 
phenomenon  of  Doppler  spread  in  which  the  spectrum  of  the  transmitted  signal  is  not 
displaced  in  the  frequency  by  Doppler  shift,  but  is  actually  spread  out  [24]. 

Doppler  spread,  is  a  measure  of  spectral  broadening,  fd  caused  by  the  time  rate 
of  change  of  the  wideband  wireless  channel.  The  Doppler  spread  is  defined  as  the  range 
of  frequencies  over  the  received  Doppler  power  spectrum  S(v).  The  components  range 
from  ifc  -fd)  to  (fc  +  fd),  where /c  is  the  transmitted  pure  sinusoidal  tone  of  frequency  and 
fd  is  the  Doppler  shift. 

The  S(v)  at  the  CR  receiver  is  then  given  by 
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S(y)  = 


f 

"  \- 

2^/to  L' 

Wm  /  J 

u. 

Otherwise 

(43) 


Where  v  is  the  Doppler  frequency  shift  and  fm  is  the  maximum  Doppler  shift 
given  by 

=  max  awD 

Where  can  be  either  positive  or  negative  depending  on  the  arrival  angle  0.  y  is 
the  relative  velocity,  and  A  is  the  signal  wavelength.  Thus,  in  this  case,  the  Doppler 
spread  is  obtained  by 

=  4  (45) 

The  knowledge  of  the  Doppler  power  spectrum  in  equation  43  allows  the 
determination  of  how  much  spectral  broadening  is  imposed  on  the  signal  as  a  function  of 
the  change  in  the  multipath  channel  state. 

Effect  of  Doppler  Shift  on  CS  algorithm  Based  Spectrum  Detection 
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Figure  48  Frequency  Reponses  PSD  Passed  Through  Fading  Channel  at  150  Max. 

Doppler  Shift 
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PSD  PSD 


Figure  49  Root  Mean  Square  Error  of  CS  Algorithms  in  (fading  plus  AWGN) 
Environment  at  150  Max.  Doppler  Shift 
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Figure  50  Frequency  Response  of  PSD  Passed  through  Fading  Channel  at  100  Max. 

Doppler  Shift 
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Figure  51  Root  Mean  Square  Error  of  CS  Algorithms  in  (fading  plus  AWGN) 
Environment  at  100  Max.  Doppler  Shift 
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Figure  52  Frequency  Response  of  PSD  Passed  through  Fading  Channel  at  10  Max. 

Doppler  Shift 
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Figure  53  Root  Mean  Square  Error  of  CS  Algorithms  in  (Fading  plus  AWGN) 
Environment  at  10  Max.  Doppler  Shift 

Figures  48,  50  and  52  consider  the  effect  of  Doppler  shift  on  the  system 
performance  at  different  maximum  Doppler  shifts.  As  the  Doppler  shift  increases,  the 
channel  varies  more  rapidly,  and  thus  the  system  performance  becomes  worse.  One  can 
also  see  that  in  the  low  SNR  region,  the  system  performance  changes  only  slightly  when 
using  the  sine  interpolator  which  is  not  related  to  channel  state  parameters.  The  results 
obtained  shows  that  though  LASSO  perform  worst  at  150  and  100  Maximum  Doppler 
shift,  the  result  at  10  maximum  Doppler  shift  for  all  the  three  algorithms  is  relatively 
close. 

However,  in  the  high  SNR  region  (SNR  >10  dB)  as  shown  in  Figures  49,  51  and 
53),  the  performance  changes  significantly.  The  reason  for  this  is  for  small  SNR  values, 
the  AWGN  is  the  dominant  factor  affecting  the  accuracy  of  channel  estimation.  However, 
in  the  high  SNR  region,  mismatch  of  the  sine  interpolator  with  the  channel  variations 
become  the  dominant  cause  limiting  the  system  performance. 

It  has  been  shown  that  it  is  possible  to  determine  the  available  spectrum  bandwidth 
using  compressive  sampling  based  signal  reconstruction  methods.  This  approach  is  robust 
whether  in  an  additive  white  Gaussian  noisy  channel  or  a  Rayleigh  fading  channel. 
Although  significantly  error  occurs  in  the  magnitude  of  the  characterization  based  on 
Rayleigh  fading  channel  plus  AWGN  in  comparison  to  just  AWGN  channel,  this  is  not 
necessarily  a  problem  since  the  goal  is  to  determine  the  occupied/unoccupied  frequency 
band  in  the  spectrum.  From  the  results  obtained  we  noticed  that  the  compression  ratio  has 
a  significant  effective  on  the  efficacy  of  CS  algorithms  used  for  dynamic  spectrum 
access.  Edge  detection  algorithms  can  be  employed  in  determining  the  unoccupied 
frequency  bands. 
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3.  ERGODIC  CAPACITY  OF  COOPERATIVE  AIRBORNE  NETWORKS  WITH 
ADAPTIVE  SOURCE  TRANSMISSION 


3,1  Methods,  Assumption  and  Procedures  for  Efficient  Computation  of  Ergodic 
Capacity  for  Cooperative  Relay  Networks 

While  the  range,  eapacity,  and  reliability  of  a  radio  link  ean  always  be  improved 
by  inereasing  transmit  power,  this  solution  is  limited  by  several  praetieal  eonsiderations: 
battery-life,  power  amplifier  size/weight,  interferenee  to  eo-ehannel  and  neighboring 
nodes,  and  probability-of-intereept  requirements,  all  render  this  strategy  ineffeetive.  An 
important  modus  operandi  to  mitigate  the  deleterious  effects  of  wireless  channel 
impairments  is  the  use  of  diversity.  Although  multiple-input-multiple  output  (MIMO) 
architectures  could  drastically  improve  the  range  and  reliability  of  beyond  line-of-sight 
and/or  over-the-horizon  communication  links  without  increasing  transmit-power  -  the 
cost,  weight,  and  poor  aero-dynamics  of  antenna  arrays  may  prohibit  their  use  on  both 
unmanned  vehicles  and  dismounted  tactical  war-fighters. 

In  multi-vehicle  cooperative  operations,  networked  nodes  in  a  tight  cluster  may 
coordinate  both  their  transmissions  and/or  receptions  to  mimic  a  space-time  processing 
system  as  if  they  were  part  of  a  single  antenna  array  platform  (e.g.,  [37]-[38]).  Thus 
cooperative  airborne  networking  could  significantly  increase  the  range  and  reliability  of  a 
long-haul  inter-cluster  communication,  thereby  improving  platform  endurance  with 
enhanced  LPI/LPD  capability,  without  using  an  antenna  array.  Several  standardization 
groups  such  as  IEEE  802.16  and  IEEE  802.11  have  also  incorporated  cooperative 
relaying  into  their  emerging  wireless  standards  (e.g..  Mobile  Multihop  Relaying  Group 
has  defined  a  multihop  relay  architecture  in  the  baseline  IEEE  802. 1 6j  standard). 

Although  an  intermediate  (relay)  node  in  a  cooperative  airborne  network  (CAN) 
may  either  amplify  what  it  receives  (in  case  of  amplify-and-forward  relaying)  or  digitally 
decodes,  and  re-encodes  the  source  message  (in  case  of  decode-and-forward  relaying) 
before  re-transmitting  it  to  the  destination  node,  we  shall  focus  on  the  amplify-and- 
forward  (analog)  relaying  scheme  in  this  report  because  it  does  not  require 
“sophisticated”  transceivers  at  the  relays.  Nevertheless,  our  framework  is  also  applicable 
for  digital  relaying  once  the  moment  generating  function  (MGE)  of  signal-to-noise  ratio 
(SNR)  is  found.  While  this  protocol  can  achieve  full  diversity  using  a  virtual  antenna 
array,  there  is  a  loss  of  spectral  efficiency  due  to  its  inherent  half-duplex  operation.  But 
this  penalty  could  be  “recovered”  to  some  extent  by  combining  the  cooperative  diversity 
with  a  link  adaptation  mechanism  wherein  the  power  level,  signal  constellation  size, 
coding  rate  or  other  transmission  parameters  are  adapted  autonomously  in  response  to 
fluctuations  in  the  channel  conditions.  Moreover,  fixed  transmission  methods  that  are 
designed  to  provide  the  required  quality  of  service  in  the  “worst-case”  scenario  are  very 
inefficient  when  better  channel  conditions  prevail. 

But  the  art  of  adaptive  link  layer  in  a  CAN  is  still  in  its  infancy  especially  when 
optimized  in  a  cross-layer  design  paradigm.  Majority  of  the  literature  on  cooperative 
diversity  were  limited  to  both  fixed  signaling  rate  and/or  constant  transmit  power  for  all 
communication  nodes.  It  is  also  not  very  clear  as  to  how  link  adaptation  (e.g.,  distributed 
power  control)  could  be  performed  with  limited  channel  side-information  (CSI),  and  what 
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are  the  benefits  of  jointly  optimizing  the  upper  layer  protoeols  with  an  adaptable  PHY 
layer?  For  instanee,  while  [39]-[41]  have  studied  the  problem  of  optimal  power  alloeation 
in  a  eooperative  wireless  network,  souree  rate-adaptation  [42]  was  not  eonsidered,  and 
more  eritieally  its  solution  requires  the  knowledge  of  CSI  of  all  links  (i.e.,  large 
overhead).  Motivated  by  these  observations,  [43]  derived  bounds  for  the  ergodie  eapaeity 
of  a  non-regenerative  adaptive-link  eooperative  diversity  system  with  limited  CSI,  in 
whieh  the  rate  and/or  power  level  at  the  souree  node  is  adapted  aecording  to  the 
prevailing  ehannel  eondition  (i.e.,  only  the  effeetive  SNR  at  the  destination  node  needs  to 
be  feedbaek  to  the  souree  node)  while  the  relays  simply  amplify  and  forward  the  signals. 
Nevertheless,  their  analysis  and  results  were  limited  to  Rayleigh  fading.  But  in  a  CAN,  it 
is  much  more  reasonable  to  model  the  channel  gain  of  each  communication  link  as  a 
Nakagami-m  or  a  Rice  random  variable  (due  to  the  increased  likelihood  of  the  presence 
of  strong  specular  components  in  an  airborne  platform). 

Therefore,  the  main  objective  of  this  study  is  to  develop  a  unified  framework  based 
of  the  MGF  method  (i.e.,  since  the  MGF  of  the  end-to-end  SNR  is  perhaps  much  easier  to 
compute  or  may  be  readily  available  owing  to  its  widespread  use  in  the  outage  probability 
and  error  rate  analysis)  for  evaluating  the  ergodie  capacity  of  cooperative  analog  relaying 
systems  under  different  adaptive  source  transmission  policies  in  a  myriad  of  fading 
environments.  Aside  from  its  utility  in  the  assessment  of  mean  achievable  transmission 
rates  in  an  adaptive-link  CAN,  our  development  is  also  interesting  because  several 
authors  [44] -[46]  have  recently  argued  that  although  the  MGF  approach  has  been 
successfully  and  extensively  applied  for  evaluating  the  performance  of  wireless  relaying 
systems  in  terms  of  outage  probability  and  error  rates,  there  have  been  very  limited 
contributions  on  ergodie  capacity  of  fading  relay  channels  [44,  pp.  2286]  (which  may  be 
attributed  to  the  difficulty  in  evaluating  the  probability  density  function  of  the  end-to-end 
SNR  in  closed-form)  or  explicitly  highlighted  the  complexity  of  generalizing  the  MGF 
approach  for  channel  capacity  computation  [45].  In  [43],  the  authors  circumvented  this 
difficulty  by  evaluating  upper  and  lower  bounds  of  the  ergodie  capacity  instead  (since  its 
PDF  for  the  Rayleigh  fading  case  can  be  expressed  in  closed-form),  while  [44]  resorted  to 
the  method  of  moments  (also  restricted  to  Rayleigh  channels).  In  [46],  the  authors 
proposed  a  general  method  for  channel  capacity  analysis  using  an  integral  relation  known 
as  Ei-transform,  but  their  solution  requires  the  derivatives  of  the  MGF  and/or  its  auxiliary 
function,  which  can  be  very  cumbersome.  Moreover,  the  ergodie  capacity  formula  for  the 
OPRA  policy  is  an  approximation  while  the  ergodie  capacity  for  the  TCIFR  policy  was 
not  considered.  In  this  report,  we  will  present  two  unified  yet  novel  and  efficient 
approaches  for  evaluating  the  ergodie  capacity  of  cooperative  relay  networks:  (i) 
cumulative  distribution  function  (CDF)  method  -  which  utilizes  the  MGF  of  end-to-end 
SNR  in  conjunction  with  a  fixed-Talbot  method  [47];  (ii)  MGF  method  -  which  exploits 
an  “exponential-type”  integral  representation  for  the  logarithm  function  that  facilitates  the 
averaging  problem  over  generalized  fading  channels. 
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3,2  System  Model 


Figure  54:  Link-adaptive  cooperative  diversity  system  model. 


Figure  54  illustrates  a  eooperative  wireless  network  model.  The  system  model  goal  is 
to  ensure  the  eonnectivity  and  network  stability  needed  to  support  varying  quality-of- 
serviee  requirements  in  an  airborne  platform.  The  souree  node  S  eommunicates  with 
destination  node  D  via  direet-link  and  through  N  amplify-and-forward  relays, 
Ri  i  e  {1,2,....,  A?},  in  two  transmission  phases.  During  the  initial  Phase  I,  S  transmits  signal 
.r  to  D  as  well  as  to  the  relays  Ri,  where  the  ehannel  fading  coeffieients  between  S  and  D, 
S  and  the  i-th  relay  node  Ri,  and  Ri  and  D  are  denoted  by  ,  a, ,  and  a.^ ,  respectively. 
The  received  signals  at  the  destination  node  and  i-th  relay  node  are  given  by 

i',,.  (46) 

y.,  i  =  .  (47) 


where  Pj  =  P ,  is  the  transmitted  power  at  the  source,  while  j ,  77,  ,.  are  the  additive 

noises  introduced  between  source  to  destination,  and  source  to  relay,  respectively.  During 
the  second  phase  of  cooperation,  each  of  the  N  relays  transmits  the  received  signal  after 
amplification  via  orthogonal  transmissions  (e.g.,  TDMA  in  a  round-robin  fashion  and  /or 
FDMA).  Suppose  P  ,,  is  the  transmit  power  of  the  i-th  relay  while  transmitting  its  signal  to 

the  destination  node.  The  received  signal  at  the  destination  during  the  second  phase  of 
operation  is  given  by  [37] 


y,a  =  +1],,  = 


(48) 


+  K 


where  G,  is  the  i-th  relay  amplifier  gain,  denoted  as  [  1  ],  G/  =  j p  .  .  f  +N„  ,  is  the 

noise  introduced  between  z-th  relay  and  destination.  We  further  assume  that  the  total 

N 

transmitted  power  P^  =  Pt  is  fixed.  Substituting  (47)  into  (48),  we  obtain 


yt.^  = 


+  N^ 


(49) 
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where 


+?;_^,is  a  zero-mean  complex  Gaussian  random  variable  with 


\a.\  +N„ 


variance 


variables. 


/ 

' 

2 

a^. 

2 

+  N^ 

+  1 


by  assuming  that  ,  and  ^  are  independent  random 


Suppose  a  maximum-ratio  combining  (MRC)  is  employed  at  the  destination  node  to 
combine  all  the  signals  received  during  Phase  I  and  Phase  II,  the  output  of  the  MRC 
detector  is  given  by 

N 


P  P. 


where  b,  = 


as  [1] 


N 


+  N 


Thus  the  total  received  SNR  can  be  shown 


Pu 

a,-,. 

2 

p. 

a^. 

2 

+  N^ 

+  1 


Yt  =  7.4  +  S 


Y 

•  s,t*  i,c 


^  +  Y.,+Yu 


i\ 

=  [  Ytb=Y.4+Yj^‘  ] 


(51) 


where  PjT  / n„  and  /v„  denote  the  instantaneous  SNR  of  the 

source-relay  and  relay-destination  links,  respectively.  The  approximation  on  the  right- 
side  of  (51)  is  obtained  by  recognizing  that  the  instantaneous  SNR  of  a  two  hops  path  can 
be  accurately  estimated  to  be  the  harmonic  mean  of  individual  link  SNRs  especially  at 
moderate  and  high  SNR  regimes  [50] -[51].  Assuming  y,  ,,  y.^  are  independent  random 

variables,  it  is  rather  straight  forward  to  show  that  the  MGF  of  y^^  can  be  evaluated  as 

('^)  =  ^r.4  ('^)n  ('^)  (^2) 

i  =  1 

Note  that  the  MGF  of  y,  =  y,.,7,4 /(y,., +f,,d)  iii  ^  Rayleigh  channel  with  independent  but 
non-identically  distributed  (i.n.d)  fading  statistics  well-known,  and  is  given  by  [50] 


(j)^^  (i)  =  [(1  /  .  - 1  /  n,  „  )N  (1  /  ,  + 1  /  n,, ,,  )i]  /  aN 


25 


( 


A’n  Q 


-In 


.  1  1 

j-i-Ah - V  — 

n  n  ,  , 

yv  "J  J 


n  n 


(53) 


where  n„,  =  £:[y„^]  corresponds  to  the  mean  link  SNR  and 
A  =  ^(1  /  n,  -  - 1  /  n,,  „  )N  25(1  /  -  -1 1  /  Q,  j  -1  5N  Moreover,  the  MGF  of  y,  in  a  Nakagami  channel 

with  independent  identically  distributed  (i.i.d)  fading  statistics  is  given  by  [51] 


^  1  Qy 

m,  2m;  m  +  —  ; - 

2  4m 


(54) 


If  a  closed-form  expression  for  y  is  not  available  (e.g.,  Nakagami -m  or  Rice  channels 
with  i.n.d  fading  statistics),  but  does  exist  for  y,=  min  (y,,,y,rf) ,  we  may  then  resort  to  the 
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development  of  eapacity  bounds  using  the  inequality  used  in  [43], 

[Tlb  =  r.,,,  +  -  X  +  X  ^  ’ 


/V 


^^{s)<c^.Xs)<c^^Js)Y\^,,{sI2) 


(i.e., 


(55) 


In  this  case,  the  MGF  of  y  for  a  2-hops  path  may  be  derived  as 


k  €  {{s,iUi,d)) 
j*k 


(56) 


For  instance,  it  is  not  difficult  to  show  that  the  MGF  of  f.  in  a  Nakagami-m  channel  with 
i.n.d  fading  statistics  is  given  by 


h  (^)  = 


T{m^+mj) 

j*k 


z 


(57) 


2-^1 


\-mj,m^;  1  +  m^; 


skliklf.  +Q^m^  +D.i^mj 


Once  the  MGF  of  y.  or  y  is  found,  we  can  compute  the  outage  probability  (i.e.,  it’s 
CDF)  efficiently  using  a  fixed-Talbot  (multi-precision  Laplace  transform  inversion) 
method  [47],  viz.. 


FXx)  = 


M 


z: 


>  ^  5(60  2M 


(58) 


where  r  =  2M /{5x),  cr(60  =  6( +(6(cot(60-l)cot(60,  =  knIM,  and 

=  ^<9,0'  +  cot(^0). 

In  Section  3.3,  it  will  be  shown  that  the  ergodic  capacities  of  ORA,  TCIFR,  and  OPRA 
source  adaptive  transmission  policies  can  be  expressed  in  terms  of  the  complementary 
CDF  and/or  the  MGF  of  end-to-end  SNR  alone.  Hence  the  upper  limit  on  bandwidth- 
normalized  link  spectral  efficiency  of  cooperative  relay  networks  may  be  evaluated 
readily  for  all  of  the  cases  discussed  above  owing  to  the  availability  of  the  MGF  of  y  or 
yin  closed-form.  For  example,  (56)  (or  (61)),  (63)  and  (66)  (or  (67))  in  conjunction  with 
(55),  (57)  and  (58)  generalize  the  results  in  [43]  to  i.n.d  Nakagami-m  channels.  Similarly, 
precise  estimates  of  the  ergodic  capacities  with  different  source  adaptive  transmission 
techniques  in  i.n.d  Rayleigh  fading  [or  i.i.d  Nakagami-m  fading]  can  be  obtained  by 
substituting  (52)  and  (53)  [or  (54)  for  i.i.d  Nakagami-m  fading]  into  (60),  (63)  and  (66). 


In  Section  3.4,  it  will  be  shown  that  the  outage  probability  can  be  expressed  in  terms 
of  CDF.  Whereas,  achievable  spectral  efficiency  and  average  bit  error  rate  can  be 
expressed  in  terms  of  difference  of  two  CDF  terms;  and/or  MGF,  in  conjunction  with 
fixed-Talbot  method.  Once  the  MGF  of  y  or  y  is  available  in  closed  from.  For  example, 
(74)  in  conjunction  with  (55)  and  (57);  and  (64)  in  conjunction  with  (52),  (53)  and  (55) 
generalize  the  results  in  [52]  for  average  spectral  efficiency  of  ACR  in  i.n.d  Nakagami-m 
channels  and  i.n.d.  Rayleigh  fading  channel  respectively.  Similarly,  accurate  estimates  of 
average  spectral  efficiency  of  ADR  in  i.n.d  Nakagami-m  channels  can  be  obtained  by 
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substituting  (76)  (by  computing  CDF  terms  using  the  fixed-Talbot  method),  in  to  (75), 
once  the  MGF  is  available  in  elosed  from.  Additionally,  average  bit  error  rate  can  be 
expressed  by  substituting  (70)  in  to  (67). 

3,3  Ergodic  Capacity  Computation  in  Fading  Channels 

The  well-known  Shannon-Hartley  law  tells  us  that  there  is  an  absolute  limit  on  the 
error-free  bit  rate  R  that  ean  be  transmitted  within  a  certain  ehannel  bandwidth  5  at  a 
specified  SNR.  This  theoretieal  limit  denotes  the  ehannel  eapaeity  C.  Shannon’s  noisy 
ehannel  eoding  theorem  also  states  that,  it  is  not  possible  to  make  the  probability  of  error 
tend  to  zero  if  R  >  C  with  any  code  design.  Thus,  it  is  elear  that  the  metrie  C  plays  an 
important  role  in  the  design  or  appraisal  of  any  communications  system  (sinee  it  serves  as 
an  upper  limit  on  the  transmission  rate  for  reliable  communieations  over  a  noisy 
communieation  ehannel).  In  this  seetion,  we  will  derive  generic  expressions  for 
computing  the  ergodie  eapaeity  of  cooperative  relay  networks  with  different  adaptive 
source  transmission  teehniques. 


3,3,1  Optimal  Rate  Adaptation  with  Fixed  Transmit  Power  (ORA) 

When  only  the  rate  is  adapted  by  changing  channel  conditions,  the  ergodic  capacity  is 
given  by 

%  =  +  (59) 

B  iV  +  lln2Jo 

Utilizing  the  “exponential -type”  integral  representation  of  ln(y  +  l)  (see  Appendix  A),  we 
ean  facilitate  the  averaging  problem  in  (48)  given  that  ^,(  )is  available  in  elosed  form. 
Substituting  (A.5)  into  (48),  we  obtain  the  ergodie  eapaeity  in  terms  of  the  MGF  of 
/  alone,  viz.. 


B 


'  N  +  1  ln2 
1  1 


dx 


(60) 


N 


Alternatively,  using  integration  by  parts  (with  u  =  \n{i  +  r),dv  =  f^ir)dr  andv  = 
(48)  can  be  expressed  in  terms  of  CDF  of  y  as 


-£  fy{y)dy). 


_ 1 

B  iV  +  lln2Jo  l  +  x 


(61) 


where  A  (x)  =  l-F  (x)  =  J  f^(r)dy  is  the  complementary  CDF  of  y. 

It  is  also  interesting  to  note  that  the  new  forms  (49)  and  (50)  allow  us  to  show  that  the 
ergodic  capacity  of  a  point-to-point  link  increases  with  the  increasing  reeeiver  diversity 
order  regardless  of  the  fading  ehannel  model  or  diversity  eombining  teehnique  employed. 
Gunther  [49,  pp.  401]  suggested  that  while  this  is  intuitive,  it  not  easy  to  prove  this  trend 
mathematieally  for  the  ORA  adaptation  policy.  However,  reeognizing  that  the  CDF 

F(7)will  deerease  with  increasing  order  of  reeeiver  diversity  for  any  given  rAhe 
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numerator  term  F‘{x)  in  (61)  will  be  much  closer  to  unity  in  comparison  with  that  of  no¬ 
diversity  or  lower  diversity  order  cases.  Thus,  we  have  shown  that  the  capacity  increases 
with  diversity  order  as  expected.  Similarly,  since  is  a  monotonically  decreasing 

function  with  respect  to  its  argument  and  \ ,  the  term  in  (60) 

approaches  to  much  closer  to  unity  as  receiver  diversity  order  increases. 


3,3,2  Channel  Inversion  with  Fixed  Rate  (CIFR) 

In  CIFR  policy,  the  transmitter  adapts  its  power  to  maintain  a  constant  SNR  at  the 
receiver  and  uses  fixed-rate  modulation  and  fixed-code  designs.  This  technique  is  the 
least  complex  to  implement  given  that  reliable  channel  estimates  are  available  at  the 
transmitter.  However,  when  the  channel  experiences  deep  fades,  the  penalty  in  transmit 
power  requirement  with  the  CIFR  policy  will  be  enormous  because  channel  inversion 
needs  to  compensate  for  the  deep  fades.  To  overcome  this,  a  truncated  channel  inversion 
and  fixed  rate  policy  (TCIFR)  was  considered  in  [42]  where  the  channel  fading  is  only 
inverted  above  a  fixed  cutoff  fade  depth  The  data  transmission  is  ceased  if  y  falls 
below  7„.In  this  case,  it  is  easy  to  show  that  the  outage  probability  is  fio)  the 

ergodic  channel  capacity  is  given  by 


C 

^  TCIFR 

B 


N  +  l 


log^ 


1  +  - 


[  r'fM)dr 

V  *'’'0  y 


^;(fo) 


(62) 


Substituting  f,(r)  =  ~\  B.e{^^{-jco)e"}do)  (i.e.,  inverse  Fourier  transform  of  the 
;r 

characteristic  function)  in  (51),  we  obtain  (see  details  in  Appendix  B) 


C 

^  TCIFR 

B 


N  +  l 


log^ 


1-- 

V 


(63) 


1  poo  ^ 

where  V  =  —  Re{^  while  the  exponential  integral  Ei{-jc)  =  -\  — dt 

can  be  evaluated  in  MATLAB  using  command  line  “cosint(c)-7(-;r/2  +  sinint(c))”. 
Observe  that,  if  the  MGF  of  y  is  known  in  closed-form,  then  the  above  integral  can  be 
evaluated  efficiently  via  Gauss-Chebyshev  quadrature  method  for  a  wide-range  of  fading 
channel  models  and  diversity  combining  techniques  employed  (including  maximal-ratio 
combining  and  selection  diversity).  The  zero-outage  ergodic  capacity  (i.e.,  CIFR  policy) 
can  be  evaluated  by  setting  =  0  in  (63). 


3,3,3  Optimal  Power  and  Rate  Adaptation  (OPRA) 

In  the  OPRA  strategy,  both  the  transmission  power  and  rate  are  matched  to  the  varying 
channel  condition  through  use  of  a  multiplexed  multiple  codebook  design.  This  leads  to 
the  highest  achievable  capacity  with  CSI.  From  [42],  we  have 


B 


1  1 
A?  + 1  In  2 


(  \ 

r_ 

Jro 

fAr)dr 


(64) 
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where  is  optimal  cut-off  SNR  below  which  the  data  transmission  is  suspended. 
Substituting  (A.4)  into  (53),  and  by  re-arranging  the  integral  terms,  we  obtain 


c,,.  _  1  1 


N  + 

1  1 

N  + 


,  r  .00 

7— [  -  e-^^\  f,{r)dr-\  e  f^(r)dr 
1  ln2  Jo  xL  J»  Jj-o 

1  ln2  Jo  x'- 


dx 


(65) 


because  f  4(y)4y  =  l-F^(y„)  =  F;(y„)  and  f  e  f^(r)dr=  ^^(2x/y„,y„)  denotes  the 

J/o  •'ro 

marginal  MGF.  Since  the  desired  marginal  MGF  is  typically  not  available  in  closed- 
form,  we  circumvent  this  difficulty  by  computing  this  term  via  the  “CDF”  of  an  auxiliary 
function  with  a  fixed-Talbot  method  (details  shown  in  Appendix  C),  viz.. 


C 

^OPRA 

B 


1  1 
W  +  1  ln2 


, 

^2x^ 

+  F(y,) 

Jo  X 

vA  J 

dx 


(66) 


If  ^^(2x/7jis  available  in  the  closed  form,  then  the  terms  F(x„)  and  F.(x„)in  (66)  can  be 

computed  efficiently  using  the  fixed-Talbot  method  (see  (58)). 

Alternatively,  (64)  can  be  concisely  expressed  in  terms  of  the  CDF  of  /  (using 

integration  by  parts  with  u  =  and  v  =  -  f  (x)dx  )  (64)  as 

^  Tor 


_ ^rElMldy. 

B  A?  +  l  ln2  J'  Y 


(67) 


To  achieve  the  capacity  (66)  or  (67),  the  channel  fade  level  (i.e.,  CSI)  tracked  at  the 
receiver  must  be  conveyed  to  the  transmitter  on  the  feedback  path  for  power  and  rate 
adaptation  in  accordance  with  the  time -varying  nature  of  the  channel  fading. 
When  y  <Y„,  no  data  is  transmitted,  and  thus  the  optimal  policy  suffers  an  outage 
probability  given  by  F„,  =  PyiYo) ,  which  equals  to  the  probability  of  no  transmission.  The 
optimal  cutoff  SNR  must  satisfy 


f;  (a  )  -  a  [1  +  f  r"fr  (Y)dY]  =  0  (68) 

The  integral  term  in  (34)  can  be  evaluated  efficiently  similar  to  the  development  of 
(63)  when  the  MGF  is  available.  Furthermore,  asymptotic  analysis  of  (68)  shows  that 
when  the  mean  SNR  n^obecause  F(x)^iand  4(x)^0(i.e.,  the  effect  of 
n  ^  0  can  be  predicted  from  the  normalized  PDF  or  the  normalized  CDF  curve  when  its 
argument  x  ^  oo ).  When  n^oo,  F(x)^0  because  this  is  equivalent  to  computing  the  CDF 
when  its  argumentx^O.  It  is  also  well-known  that  ^0 asn^oo .  Hence, as 
Q  ^  00 .  Thus  (can  be  determined  by  solving  (68)  numerically)  always  lies  in  the 
interval  [0,  1]  regardless  of  the  assumption  on  the  fading  channel  model  or  the  number  of 
relay  nodes  employed. 
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3,4  Adaptive  M-QAM  Modulation 
3,4,1  Adaptive  Scheme 


It  has  already  been  proved  that  optimal  rate-fixed  power  adaptation  policy  has 
numerous  practical  advantages,  in  addition  to  other  adaptive  transmission  schemes. 
Adaptive  QAM  method  is  commonly  used  to  enhance  the  link  spectral  efficiency  of 
system,  by  satisfying  a  certain  required  BER  performance  [53]-[56];  and  in  [53]  the 
authors  nicely  propose  and  analyze  a  variable -rate  variable  power  M-QAM  for  fading 
channels.  By  considering  above  facts,  adaptive  M-QAM  is  considered  in  this  paper.  The 
equation  of  bit  error  rate  (BER)  of  a  system,  which  incorporate  M-QAM  modulation  over 
an  Additive  White  Gaussian  Noise  (AWGN)  channel,  with  coherent  detection;  and  Gray 
coding  can  be  represented  as  [53],  [56] 


BER{M^ 


Y)  *0.2  exp 


3y 


2(M-1) 


(69) 


where  yis  the  received  SNR.  As  from  [53],  [56],  (69)  provide  closed-form  solution  for 
further  analysis,  and  it  is  invertible.  By  inverting  (69),  the  spectral  efficiency  of  the 
continuous-rate  M-QAM  can  be  approximated  as  [53],  [56]; 


R 

IT 


=log,M=log, 


1  + 


3y 


2K, 


(70) 


0  J 


where  =  -^(Sfi^R,)  and  BER^  represents  the  target  BER. 

We  can  effectively  compute  the  spectral  efficiently  of  Adaptive  Continuous  Rate  M- 
QAM  using  (59).  However,  compared  to  ACR  M-QAM,  ADR  M-QAM;  in  which  the 
constellation  size  is  limited  to  2"  (where  n  is  positive  integer),  is  more  realistic  approach. 

In  ADR  M-QAM,  The  range  of  the  effective  received  SNR  is  divided  into  T+1  fading 
regions.  In  each  region,  a  specific  constellation  of  size  is  used.  A  particular 
constellation,  of  size  is  decided  at  the  receiver,  based  on  the  overall  total  received 
SNR,  when  the  fading  causes  the  effective  SNR  to  fall  into  the  nth  region  (n  =  0,  1,  .  .  .  , 
T).  The  partitioning  of  the  effective  received  SNR  depends  on  the  desired  BER 
level and  to  reach  a  specific  BER  target  of  BER^^,hy  using  m„-QAM;  region 
boundaries  are  selected  to  the  following  range  of  SNR,  to  accomplish  desired 
performance  [56]; 

y  =[erfc-'{2BERJf, 


y„  =-ir„(2"-l);  «  =  0,2,3,....,r, 


(71) 


A+i=+“ 

where  e//c  '  ()  is  the  inverse  complementary  error  function. 


B.  Probability  of  Outage 

Wheny<y, ,  no  data  is  transmitted,  and  thus  the  optimal  policy  suffers  an  outage 
probability,  given  by  f,,,,  =  f,(yi)  where. 
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F,{r,)  =  LT- 


(72) 


where  7,  is  the  optimal  eutoff  SNR  level,  below  which  data  transmission  is  suspended. 
Once  the  MGF  of  y.  or  is  found,  we  can  compute  the  outage  probability  (i.e.,  it’s 
CDF)  efficiently  using  a  fixed-Talbot  (multi-precision  Laplace  transform  inversion) 
method  [47]  (see  (47)). 


C.  Achievable  Spectral  Efficiency 

The  average  spectral  efficiency  for  ACR  is  calculated  by  integrating  (59)  over  the  PDF 
of  the  received  SNR  [52]. 


Kcr 

B 


_[ _ 1_ 

N  +  \  In  2  Jo 


In 


^+j^\frir)dr 


(73) 


3r_ 

2K„ 


(see  Appendix  A), 


Utilizing  the  “exponential  type”  integral  representation  “of  In  1  + 

I 

we  can  efficiently  compute  achievable  spectral  efficiency  given  that,  ^^(Ois  available  in 
closed  form.  Substituting  (A.5)  into  (52),  we  obtain  the  average  spectral  efficiency  in 
terms  of  MGF  of  y  alone,  viz.. 


1  1 


Fgc,  ^ _ 

B  N  + 

1  1 


N- 


N  + 


1  In  2  Jo  X  Jo  ' 

1  r”  \  (15) 

TTA  — 

lln2Jo  X  l-^oJ 

dy 


dx 


dx 


(74) 


1  1 


_^r 

1  In  2  J » 


In  comparison  to  achievable  spectral  efficiency  of  ACR,  the  achievable  spectral 
efficiency  of  ADR  is  defined  as  the  sum  of  data  rates  in  each  of  the  partition  regions, 
weighted  by  the  probability  of  occurrence  of  each  region,  and  is  given  by[52]. 


R 


B 


(75) 


N  +  lf-f 

where  n  =  log2(MJ  corresponds  to  the  data  rate  of  the  nth  region,  is  denoted  as 
the  probability  of  occurrence  that  the  effective  received  SNR,  y  is  in  the  nth  partition 
region  and  is  given  as[16]. 


fii+1 

p„  =  j  f  (y)dy  =  F  (y„„ )  -  f  (y„ ) 

rn 


(76) 


Here,  the  MGF  of  random  variable  y  is  available  and  by  (47),  Abate ’s  fixed  Talbot 
method  (multi-precision  Laplace  inversion  formula)  [47],  we  can  efficiently  evaluate  the 
CDF  terms  F^(y„^, )  andF^(yJ. 


D.  Average  Bit  Error  Rate 
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The  analysis  of  is  insignificant,  because  the  value  of  ACR  M-QAM 
constantly  approaches  the  desired  BER  target.  However,  in  the  case  of  ADR  M-QAM,  the 
received  SNR  is  divided  in  to  discrete  partitions,  results  in  a  conservative  average  BER. 
This  limits  the  average  to  be  always  smaller  than  the  target  BER.  The  average 

BER^^^^  can  be  approximated  as  [52], 


BER 


n  BER^ 


n 


(77) 


r,i+i 

where  BER=  j  BERiM^^  r)frir)dr 

Yn 

by  substituting  (14)in  (23)  we  get, 

Yn+l  / 

J"  0.2  exp 


BER  = 


3/ 

2(M-1) 


=  0.2 


jexp 


3r 


2(M  -l) 


frir)dr 


fy{r)dr-  I 


exp 


2(M-1) 


L(r)dr 


(78) 


(79) 


where,  f.(x)=a^  exp(-;0x)/^(x)is  an  auxiliary  function,  in  which  = 


3 

2(M  -1) 


and 


=  0.2. 

By  using  (C.l)  (see  Appendix  C),  we  can  re-write  (45)  in  terms  of  marginal  MGE,  and 
subsequently  in  terms  CDE  as  follows, 

^  =  o.2{^^(Ay,)-?5,(Ar„J} 

=  0.2  |(k  (00)  -  F.  (r„ ))  -  (k  (00)  -  K  (y_^,  ))|  (80) 

=  0.2{F.ir„,,)-F.ir„)} 


Here,  the  MGE  of  new  random  variable  f  is  available  and  by  using  (C.2),  in 
conjunction  with  (58)  Abate’s  fixed  Talbot  method  (multi-precision  Laplace  inversion 
formula)  [47],  we  can  efficiently  evaluate  the  CDF  terms  E(y„+, )  andE.(yJ  . 

Finally,  using  (76)  and  (80),  in  to  (77),  we  can  efficiently  compute  the 
average  BER^^^ . 

3,5  Results  and  Discussions  of  Ergodic  Capacity  of  Cooperative  Networks 

3,5,1  Ergodic  Capacity  of  Cooperative  Relay  Networks 

In  this  sub-section,  selected  numerical  results  will  be  presented  for  the  ergodic 
capacities  of  cooperative  analog  relaying  systems  under  different  source  adaptive 
transmission  policies.  The  following  mean  link  SNRs  (arbitrarily  chosen)  will  be  used  to 
generate  the  plots,  unless  stated  otherwise:  =  E/No,  0.5E/No,  0.5E/No, 

E/Nq,  and  =  0.2E/No. 

Figure  55  depicts  the  comparison  of  ergodic  capacities  of  three  distinct  source  adaptive 
transmission  policies  in  an  i.n.d  Rayleigh  channel  with  two  cooperative  relays.  As 
anticipated,  there  is  no  significant  difference  observed  in  the  ergodic  capacities  of  OPRA 
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and  ORA  at  high  SNRs.  But  the  ergodic  capacity  of  TCIFR  policy  (plotted  for  the  cut-off 
SNR  Xo  =  6  dB)  is  considerably  lower  than  the  OPRA  scheme.  Although  not  shown  in 
this  figure,  we  also  noticed  that  the  curves  corresponding  to  the  “tight  approximation” 

N 

case  (i.e.,  )  are  in  good  agreement  with  the  Monte  Carlo  simulation  results. 

Moreover,  the  actual  ergodic  capacity  is  slightly  closer  to  the  lower  bound  (rather  than 
the  upper  bound)  at  low  SNRs.  Although  the  authors  in  [43]  have  studied  the  channel 
capacities  of  cooperative  relaying  system  in  i.n.d  Rayleigh  channel,  but  their  framework 
does  not  lend  itself  to  the  analysis  of  the  “tight  approximation”  case  or  generalize  to  other 
fading  channels,  whereas  our  framework  encapsulates  all  these  cases  in  an  unified  way. 


V  Tight-Approximation 
□  Upper-Bound  (y^g) 

o  Lower-Bound  (y|_g) 

- ORA 

-  -  -  TCIFR  (y^  =  6dB) 
- OPRA 


5  10  15 

SNR  Es/No  (dB) 


Figure  55  Ergodic  capacities  of  ORA,  OPRA  and  TCIFR  policies  in  an  i.n.d 

Rayleigh  channel  (N  =  2). 
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Cut-Off  SNR  (dB) 


Figure  56  Ergodic  capacity  of  TCIFR  versus  cut-off  SNR  (^g)  in  an  i.n.d  Rayleigh 

channel  (N=2). 

Figure  56  depicts  the  ergodic  capacities  of  TCIFR  scheme  for  E/Nq  =  6dB  and  E/Nq  = 
15dB  as  a  function  of  cut-off  SNR  (xj  in  an  i.n.d  Rayleigh  channel  with  two  relays.  It  is 
evident  that  there  exists  an  optimal  choice  for  the  cut-off  SNR  which  maximizes  the 
channel  capacity  when  E/No  is  fixed.  But  it  should  be  also  emphasized  that  the  selection 
of  will  directly  affect  the  outage  probability  (i.e.,  probability  of  no  transmission)  of  the 
TCIFR  policy. 

Figure  57  shows  the  ergodic  capacities  (i.e.,  for  the  “tight-approximation”  case)  for 
different  source  adaptive  transmission  schemes  over  i.i.d  Nakagami-m  channels  (fading 
severity  index  m  =  0.5,  1,  1.5  and  2).  It  is  observed  that  the  ergodic  capacity  increases 
with  the  increasing  value  of  m  (i.e.,  as  channel  experiences  less  severe  fading)  for  all 
source  adaptive  transmission  schemes.  However,  the  gap  between  the  curves 
corresponding  to  the  ORA  and  OPRA  policies  widens  as  m  decreases,  although  their 
achievable  normalized  average  transmission  rates  is  quite  similar  at  high  F/Ao  values. 
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Figure  57  Ergodic  capacities  of  ORA,  OPRA  and  TCIFR  policies  in  an  i.i.d 

Nakagami  channel  (N  =  1). 


Figure  58  depicts  the  ergodic  capacity  of  ORA  in  a  Nakagami-m  channel  with  one 
relay  diversity  (N  =  1)  in  comparison  to  no  relay  case  (N  =0  )  and  for  different  transmit 
power  assignments.  Without  any  loss  of  generality,  assume  that  the  mean  SNRs  on 
different  links  are  given  as  follows  (which  takes  into  account  of  the  relative  distances 
among  the  cooperating  nodes): 


(d  ^ 

h 

(d  ) 

a,= 

( A  y " 


Ej./N  ,and  =  —  c5  ^E^jN  ,  where  c  is 
constant  that  is  related  to  the  carrier  wavelength  (i.e.,  c  =  [  —  ]  for  free  space  path  loss), 


p  p 

Lp  is  the  path  loss  exponent,  5^^  =S,.  =-^  while  =(1-^,,)  =  -^  since  =1,  and 

denotes  the  distance  between  link  a  -  b.  To  generate  Figure  58,  we  have  arbitrarily 
chosenc?^j  =d^J  =500m,  c  =  10'  =  1000m  and  L^=  4.  It  is  evident  that  equal  power 

allocation  strategy  may  be  reasonable  but  not  optimal  in  all  cases  (i.e.,  optimum  power 
allocation  strategy  is  more  beneficial  when  the  channel  experience  more  severe  fading). 
We  also  observe  that  the  tightness  of  upper  and  lower  bounds  on  the  ergodic  capacity 
becomes  looser  as  m  increases.  Comparison  between  the  curves  corresponding  to  the 
single  relay  and  no  relay  cases  reveal  that  the  optimum  power  assignment  is  not  strongly 
influenced  by  the  source-destination  link  but  rather  on  the  link  qualities  of  source-relay 
and  relay-destination  links. 
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Figure  58  Ergodic  capacity  of  N  relays  with  ORA  policy  as  a  function  of  fixed 
transmit  power  assignment  ,  in  a  Nakagami-m  channel. 


Figure  59a:  Ergodic  capacity  (“Upper  Bound  Case”)  of  N  relays  (N=0, 1)  with  ORA 

policy  in  a  Nakagami-m  channel 
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Figure  59b:  Different  transmit  power  assignment,  S^  .  (“Upper  Bound  Case”)  of  N 
relays  (N=0, 1)  with  ORA  policy  (c  =  10'^) 


Figure  59c:  Ergodic  capacity  (“Upper  Bound”  Case)  of  N  relays  (N=0, 1)  with  ORA 
policy  in  a  Nakagami-m  (m=1.5)  fading  channel. 


Figure  59a  shows  the  Ergodie  eapaeity  of  ORA  poliey  in  Nakagami-m  ehannel  as  a 

d  . 

function  of  ratio  of  distances  between  source  to  relay  and  source  to  destination  (i.e.,  — ^  ), 

d  , 

s,a 

which  describe  the  important  scenario  of  effect  of  node  placement  in  cooperative 
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communication.  For  simplicity,  explained  here  for  only  one  relay  case  but,  this  analysis 
can  be  easily  extended  for  any  number  of  relays.  The  mathematical  expressions  and  other 
parameters  are  assumed  to  be  same  as  explained  above  for  Figure  58,  except  c  and  the 
distances.  We  have  arbitrarily  chosen  values  of  d^.  varied  in  the  range  of  0  -  1000  m  with 

step  size  100  m,  c  =  lO'^d^^  =  1000m,  and  1^=  4.  It  is  observed  from  the  curves  that,  the 

position  of  relay  near  to  the  destination  (i.e.,  for  the  distance  ratio  above  0.5)  gives 
slightly  higher  spectral  efficiency  compared  to  the  values  when  the  relay  near  to  the 
source  (i.e.,  for  the  distance  ratio  below  0.5).  Furthermore,  comparison  between  curves 
corresponding  to  the  no  relay  and  single  relay  cases  disclose  that,  for  lower  and  medium 
values  of  fix  transmitted  power  (i.e.,  Et/Nq  =  5  dB  and  15  dB  respectively),  the  role  of 
relay  is  very  important  for  chosen  value  of  distance  ratio,  to  achieve  the  higher  spectral 
efficiency.  But  for  the  higher  values  of  Ej/Nq  (i.e.,  25  dB  and  above)  it  is  obvious  that  the 
use  of  relay  is  not  worth. 

Figure  59b  shows  the  Different  transmit  power  assignment,  5^.  (for  “Upper  Bound 
Case”)  of  a  cooperative  relay  network  consisting  of  N  relays  (N=0,  1)  with  ORA  policy 


as  a  function  of  ratio  of  distances,  — ^  in  a  Nakagami-m  channel  with  same  parameters  as 

d  , 

s,d 

mentioned  for  Figure  59a.  Furthermore,  comparison  between  curves  corresponding  to  the 
no  relay  case  and  single  relay  cases  reveal  that,  with  the  single  relay  case,  the  distribution 
of  transmit  power  is  distance  based.  Which,  indicates  that,  single  relay  case  utilize 
optimum  power  compared  to  the  no  relay  case. 


Figure  59c  describe  the  related  scenario  as  explained  in  Figure  59a  but,  for  the  different 
values  of  c  with  Nakagami  fading  channel  (m=1.5).  It  is  observed  from  the  curve  that,  as 
the  distance  between  the  source  and  destination  increases  (i.e.,  c  increases  from  10"'  to  10' 
)  relay  plays  an  important  role  to  achieve  desirable  spectral  efficiency.  Finally  From 
Figure(s)  59a,  b  and  c,  we  can  firmly  conclude  that,  the  spectral  efficiency  is  not  only  a 
function  of  (i)  Ej/Nq,  (ii)  fix  transmit  power  assignment  5^  ■ ;  and  (iii)  fading  index  m,  but 

also  affected  by  the  position  of  relay  (node  placement)  with  respect  to  source  and 
destination  and  distance  between  source  and  destination  (i.e.,  function  of  c). 
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1.4 


Figure  60  Ergodic  capacities  (“Upper  Bound  Case”)  in  an  i.n.d  Nakagami-m  fading 

channel  with  N  relays  (N=0, 1) 


Figure  60  shows  the  comparison  of  ergodic  capacities  of  three  distinct  source  adaptive 
transmission  policies  in  an  i.n.d  Nakagami-m  fading  channel  consisting  of  N  relays  (N=0, 
1)  with  the  following  values  of  fading  index  on  the  corresponding  link 


denotes  the  fading  index  between  link  a -b.  Without  any  loss  of  generality,  assume 

that  the  mean  SNRs  on  different  links  are  given  as  follows  (which  takes  into  account  of 
the  relative  distances  among  the  cooperating  nodes); 


Q.,  = 


fd  ^ 

-h 

(d  ^ 

h 

fd  ^ 

"..,1 

cry„.£p/fV„,Q.,,= 

[d.J 

cS,^Ej.lN^,  and 

^s,d 

constant  that  is  related  to  the  carrier  wavelength  (i.e., 


1 


where  c  is 
for  free  space  path  loss), 


P  P 

Lp  is  the  path  loss  exponent,  =  J,,  =-^  while  S.,  =(1-^,,)  =  —  since  =1,  and 


j  denotes  the  distance  between  link  a  -  b.  To  generate  Figure  60,  we  have  arbitrarily 
chosen I  =  400m,  =  600m,  c=10'^,  =  1000m  and  Lp=4.To  plot  for  N  =  1,  the 

parameters  chosen  are  <5^  j  =0.28,  <5,^  =  0.72,  <5^  =0.28  for  optimal  power  allocation 
strategy  ( )  and<7  ,  =0.5,  =  0.5,  <7  ^  =0.5 ,  for  equal  power  allocation  strategy 
i^equai )'  ^^r  thc  curvcs  of  N=  0,  the  value  of  =1  for  ORA,  OPRA  and  TCIFR. 
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SNR  E-^Nq  (dB) 

Figure  61a  Outage  probability  (“Upper  Bound  Case”)  of  source  adaptive  policies 
in  an  i.n.d  Nakagami-m  fading  channel  with  N  relays  (N=0, 1). 

Figure  61a  shows  the  eomparison  of  probability  of  outage  for  the  OPRA  and  TCIFR 
policies.  Figure  61a  is  plotted  based  on  the  cut-off  SNR  values  from  Figure  61b  for 
OPRA  and  TCIFR  as  a  function  of  SNR.  In  the  case  of  TCIFR,  the  optimal  cut-off  SNR 
is  selected,  which  maximize  the  capacity  at  the  cost  of  increased  probability  of  outage. 
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Figure  61b  Cut-off  SNR  ( 7„)  versus  SNR  Es/No  (“Upper  Bound  Case”)  in  an  i.n.d 

Nakagami-m  fading  channel 


3,5,2  Adaptive  M-QAM  Modulation  in  Cooperative  Relay  Neworks 

In  this  sub-section,  we  have  presented  numerical  results  for  the  outage  probability, 
achievable  spectral  efficiency,  and  BER  analysis  for  cooperative  systems  with  adaptive 
M-QAM  transmissions  under  Nakagami-m  channels.  For  all  the  numerical  results,  the 
system  model  as  in  [52]  is  considered.  However,  the  authors  in  [52]  have  presented  the 
above  numerical  results  of  cooperative  relaying  in  i.i.d  Rayleigh  fading,  but  their  frame 
work  does  not  lend  itself  to  the  analysis  of  the  “tight-approximation”  case  or  generalized 
to  other  fading  channels(i.e.,  Nakagami-m  or  Rice),whereas  our  framework  encapsulates 
all  these  cases  in  an  unified  way.  In  our  analysis  all  the  numerical  results  for  the  case  of 
non-i.i.d.  Nakagami-m  fading,  are  obtained  with  two  relays.  The  averages  SNR  of  the 
links  are  chosen  arbitrarily  as  in  [52]  such  that  it  represents  a  realistic  model  of  a  practical 
cooperative  communication  system.  They  are  as  follows:  Q^^=E/No,  O.8E/N0, 

Q,^=  Q.3E/No,  0.56Es/No,  and  0.2Es/No. 

In  all  of  the  non  i.i.d  cases,  the  fading  index  over  each  link  is  set  as  per  the  following  3 
Cases  with  2  relay  diversity  only. 

Case  l:m^j=l,m^2  =  =1,  Case  2:  m,,  =4, =  =  =1; 

Case  3:  =1.  where  denotes  the  fading  index  between 

link  a-b 
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Figure  62  Outage  probability  of  i.n.d.  Nakagami-m  fading  channels  with  upper 

bound  for  2  relay  diversity. 


Figure  63  Outage  probability  of  i.i.d.  Nakagami-m  fading  channels  for  upper  bound 

for  2  relay  diversity. 
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Figure  62  show  the  probability  of  outage  eomputed  using  the  numerieal  method  for 
inversion  of  Laplaee  transform  (17).  To  validate  the  analysis,  Case  1  exaetly  matches  the 
results  in  [52]  for  non-i.i.d  Rayleigh  fading  case.  From  figure,  it  can  be  concluded  that, 
for  a  given  different  average  SNRs  over  each  relay  link,  links  with  better  channel 
conditions  (as  a  result  of  higher  Nakagami  fading  index)  performs  better  than  the  case; 
where  all  the  links  have  the  same  fading  index  with  m=l.  It  is  also  observed  that  outage 
probability  performance  of  case  3  is  better  than  that  of  case  2.  Case  3  has  higher  fading 
index  in  the  relay-to-destination  links,  compared  to  source-to-destination  links.  The 
average  SNR  on  the  relay-to-destination  link  are  lower  than  the  source-to-relay  link.  This 
suggests  that  a  cooperative  communication  system  can  perform  better  with  links,  that 
experience  better  fading  channels,  even  under  lower  average  received  SNRs.  This  means 
that  transmitter  power  allocation  over  the  relays  can  be  adjusted  to  reduce  power 
consumption  by  knowing  the  channel  quality.  If  the  fading  index  in  the  Nakagami-m 
channel  is  known  to  be  higher  than  the  power  allocated  to  the  transmitter,  in  those  links  it 
may  be  reduced  without  affecting  the  overall  performance  of  the  system. 

Figure  63  depicts  the  outage  probability  of  i.i.d.  Nakagami-m  fading  channels  for 
upper  bound  along  with  special  case:  Rayleigh  fading  channel  (i.e.,  m=l),  there  is  a 
surprising  effect  of  increasing  the  Nakagami-m  fading  index  on  all  the  links.  The  average 
SNRs  over  each  link  is  kept  the  same.  The  fading  index  over  all  the  links  are  also  kept  the 
same.  The  diversity  is  provided  with  2  relays.  It  can  be  observed  that  with  increasing 
fading  index,  the  outage  probability  performance  improves  drastically  at  higher  SNRs 
(SNR  >  10).  The  above  analysis  also  concludes  that,  our  analysis,  is  viable  for  the  non¬ 
integer  fading  index  Nakagami-m  (i.e.,  m=0.5,  1.5  or  2.5)  values. 


Figure  64  Average  BER  for  i.n.d.  Nakagami-m  fading  channel  for  upper  bound  for 

2  relay  diversity. 
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Figure  65  Average  BER  for  i.n.d.  Rayleigh  fading  channel  for  2  relay  diversity. 

Figure  64  shows  the  BER  eomputed  using  (23).  The  performanee  under  ease  2  and  3 
are  better  than  ease  1.  Performance  under  case  3  is  slightly  better  than  case  2.  This  is 
because  a  higher  fading  index  in  links  with  lower  average  SNRs  improves  end-to-end 
reliability.  For  links  with  higher  SNR  improvement  in  channel  condition  does  not  affect 
reliability  significantly. 


Approved  for  Public  Release;  Distribution  Unlimited. 


77 


Figure  66  Achievable  rates  for  i.n.d.  Nakagami-m  fading  channels  for  2  relay 

diversity. 


Figure  67  Achievable  rates  for  i.n.d.  Rayleigh  fading  channels  for  2  relay 

diversity. 
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Figure  66  shows  the  achievable  spectral  efficiency  for  adaptive  continuous  rate  M- 
QAM  (19)  and  adaptive  discrete  rate  M-QAM  (20)  systems  for  target  BER  level  of 
BER^^  =10“ffior  the  non  i.i.d  fading  model.  The  plot  shows  the  comparison  of  three  cases 
of  Nakagami-m  fading  channels,  along  with  the  special  case  m=l  (Rayleigh  fading). It  can 
be  observed  that  the  capacity  per  unit  bandwidth  for  case  2  and  case  3  are  higher  than 
case  1.  The  difference  in  capacities  between  case  2  and  3  are  negligible.  Hence  channel 
conditions  do  not  have  much  effect  on  capacity.  Adaptive  discrete  rate  M-QAM  suffers  at 
least  an  additional  1 .5  dB  compared  to  the  spectral  efficiency  of  Adaptive  continuous  rate 
M-QAM. 


4.  CONCLUSIONS 

In  section  2.2,  our  investigation  showed  that  although  CS  measurements  may  be 
performed  with  no  apriori  knowledge  of  the  signal  (other  than  the  fact  that  the  signal  is 
sparse  or  compressible  in  some  basis),  the  reconstruction  error  can  depend  strongly  on  the 
choice  of  basis  function.  Hence  a  means  of  determining  the  best  sparsifying  basis 
function  is  necessary.  We  were  able  to  demonstrate  that  it  is  possible  to  predict  the  best 
sparsifying  bases  for  a  compressive  sensing  based  signal  reconstruction.  Our 
investigation  revealed  that  the  Gini  index  and  the  /o-norm  measure  are  a  good  predictor  of 
best  sparsifying  bases.  In  section  2.3,  we  demonstrated  that  the  CS  based  NYFR  is 
capable  of  detecting  signals  in  the  spectrum  without  apriori  information.  The  NYFR 
being  a  CS  based  system  requires  sub-Nyquist  sampling  rate  to  detect  available  signals  in 
the  spectrum.  The  accuracy  of  the  prediction  is  very  high  and  the  automatic  detection 
algorithm  which  was  developed  eliminates  error  associated  with  visual  analyzes  of  the 
result  and  can  aid  in  automated  and  embedded  system.  The  system  is  robust  both  in 
additive  Gaussian  noise  as  well  as  in  fading  environment.  Results  obtained  from  our 
discussion  in  section  2.4  reveals  that  it  is  possible  to  determine  the  available  spectrum 
bandwidth  through  Compressive  Sensing  based  signal  reconstruction.  This  is  achievable 
with  or  without  noisy  channel  condition.  With  the  help  of  the  results  of  section  2.2,  we 
are  able  to  determine  the  best  sparsifying  basis  for  achieving  good  reconstruction  of  the 
Power  spectrum  density  is  the  Haar  wavelet.  The  PSD  is  analogous  to  a  block  type  signal 
and  as  such  the  Haar  wavelet  is  the  best  sparsifying  basis  for  use.  Compressive  sensing 
based  signal  reconstruction  is  scalable  in  its  application.  One  can  quickly  scan  through  a 
broad  frequency  spectrum  to  determine  availability  and  subsequently  narrow  down  on  a 
section  of  interest.  We  have  investigated  the  effect  of  sample  size,  compressive  sensing 
algorithm;  on  the  quality  of  results  as  well  as  the  effect  of  different  channel  noise.  Our 
results  reveal  that  OMP  performs  better  with  higher  sampling  size  as  compared  to  BP  and 
FASSO.  LASSO  algorithm  performs  better  as  the  sample  size  is  reducing.  However  since 
the  goal  is  to  determine  the  occupied  frequency  band,  our  investigation  revealed  that  the 
CS  based  method  can  efficiently  reconstruct  the  PSD  for  adequate  detection  of  occupied 
frequency  band. 

In  section  3,  we  have  derived  the  ergodic  capacity  bounds  for  non-regenerative 
(amplify-and-forward)  cooperative  relay  networks  with  adaptive  source  transmission. 
Unlike  previous  related  studies,  we  advocate  a  simple  numerical  procedure  for  unified 
analysis  of  ergodic  capacity  of  source  adaptive  relay  networks  in  a  myriad  of  fading 
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environments.  Specifically,  two  unified  approaches  have  been  presented  for  the  ergodic 
capacity  analysis  over  generalized  fading  channels:  (i)  moment  generating  function 
method;  (ii)  cumulative  distribution  function  method.  This  framework  allows  us  to  gain 
insights  as  to  how  fade  distributions  and  dissimilar  fading  statistics  across  the  diversity 
paths  affect  the  maximum  mean  achievable  transmission  rate,  without  imposing  any 
restrictions  on  the  fading  parameters.  Even  for  the  specific  case  of  Rayleigh  fading,  our 
framework  lends  itself  to  the  development  of  a  much  tighter  bound  for  the  ergodic 
capacity  of  source  adaptive  relay  networks.  Moreover,  our  analytical  framework  also 
facilitates  the  investigation  of  optimal  power  assignment  and  node  placement  (i.e., 
position  of  relay  with  respect  to  source  and  destination)  in  cooperative  relay  networks. 
The  three  adaptive  source  transmission  techniques  considered  were:  (i)  optimal 
simultaneous  power  and  rate  adaptation  (OPRA);  (ii)  constant  power  with  optimal  rate 
adaptation  (ORA);  (iii)  truncated  channel  inversion  with  fixed  rate  (TCIFR).  Asymptotic 
bound  is  also  derived  which  shows  that  constant  power  with  optimal  rate  adaptation 
policy  provides  roughly  the  same  ergodic  capacity  as  the  optimal  simultaneous  power  and 
rate  adaptation  policy  at  high  SNR  regime.  The  gap  between  ORA  and  OPRA  curves  also 
becomes  negligible  as  the  channel  experience  less  severe  fading  (higher  fading  index  m). 
The  new  framework  may  be  employed  to  solve  a  number  of  open  research  problems  in 
adaptive  relay  networks  such  as  node  selection  strategies,  ergodic  capacity  analysis  of 
decoding-and-forward  relaying,  resource  allocation  schemes  among  cooperating  nodes, 
and  the  efficacy  of  distributed  space-time  coding. 
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APPENDIX  A 

In  this  appendix,  a  sketch  of  our  derivation  for  an  “exponential-type”  integral 
representation  of  Iny  is  provided.  Such  a  representation  will  facilitate  the  averaging 
problem  that  typically  encountered  in  capacity  analysis  over  fading  channels,  and 
therefore  leads  to  a  unified  approach  for  calculating  the  ergodic  capacity  of  CANs  in  a 
myriad  of  fading  environments  and  relaying  strategies. 

Utilizing  [12,  eq.  (1.512.2)],  we  have 
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where  y  = .  Substituting  /=—!—(’  [12,  eq.  (3.381.4)]  into  (A.  1),  we  obtain 
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^^x‘ ‘  [12,  Eq.  (1.411.2)],  (A.2)  can  be  re-stated  as 
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Finally  using  variable  substitution  x  =  z(r-l),  dz= - ,  we  arrive  at  (A.4)  after  some 
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routine  algebraic  manipulations: 


It  is  also  obvious  that 
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APPENDIX  B 

Let  e  “A(x)dxand  ‘l>x(7'^)=f  denote  the  moment  generating 

function  (MGF)  and  the  characteristic  function  (CHF)  of  random  variable  Z>0, 
respectively.  In  this  case,  the  CHF  is  related  to  the  MGF  as  O^iycj)  =  A(-y^y)  •  The  PDF 
of  X  (may  be  expressed  as  an  inverse  Fourier  transform  of  its  CHF)  is  given  by 

fx  (•^)  =  :^  [  <!>.  (-jo))e~  '‘“'‘d CO  (B .  1 ) 

Let  us  express  the  CHF  of  random  variable  X  in  its  polar  formO^);®)  =  |<l>^(;®)|e^''‘‘“\ 
Hence  (B.l)  may  be  re-stated  as 
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with  the  aid  of  identities  (B.2)  and  (D.l). 
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APPENDIX  C 

It  is  important  to  note  that  the  knowledge  of  marginal  MGF  of  end-to-end  SNR  may  be 
required  while  evaluating  the  ergodie  eapaeity  with  OPRA  poliey  (e.g.,  see  (20)). 
However,  this  quantity  is  generally  not  available  in  elosed-form.  But  if  a  elosed-form 
expression  for  the  MGF  (/)^(-)  is  available,  we  may  then  use  Abate ’s  fixed  Talbot  method 
(multi-preeision  Laplaee  inversion  formula)  for  eomputing  the  desired  marginal  MGF 
effieiently  as  highlighted  below. 


Let  us  define  an  auxiliary  funetion  f.{x)  =  e^^‘f^{x) .  Henee, 
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where  FAy)=  f  f.{x)dx=  f  e^'^"f  ix)dx.  It  is  obvious  that 
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Furthermore,  we  have 

(C.2) 

^^^^(■?)  =  1^  e^'"f-{x)dx. 

ff^^)dx, 

=  <l>yis  +  P) 

(C.3) 

Substituting  (C.2)  and  (13)  (in  eonjunetion  with  (C.3))  into  (C.I), 
numerieally  effieient  method  for  evaluating  the  required  marginal  MGF 
fading  ehannel. 

we  arrive  at  a 
in  a  generalized 
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APPENDIX  D 

While  the  exponential-integral  Ei{x)  is  usually  defined  for  real  x  <  0 ,  in  the  following  we 
will  show  that  this  funetion  is  also  well-defined  even  if  its  argument  is  purely  imaginary. 
This  is  partieularly  interesting  beeause  our  unified  expression  for  the  ergodie  eapaeity  of 
CANs  with  truneated  ehannel  inversion  poliey  and  the  transeendental  equation  for 
eomputing  the  optimal  eut-off  SNR  for  the  OPRA  poliey  are  expressed  in  terms  of 
Ei{-  jc)  where  c  >  0  is  real. 

Letting  q  =  ±jy{y  is  real),  we  have 

f”  e“'  f”  e~‘‘' 

Ei{-q)  =  -^  — dt  =  —^  - dt  (D.l) 


Utilizing  the  Euler  identity  e*"  =cos(y)  +  ysin(y),  (D.l)  ean  be  re-stated  in  terms  of  the 
more  familiar  sine-integral  and  eosine-integral,  viz., 

EK±jy)  =  -[  -cos(yOd?±f  -s\n{yt)dt, 

f  t  (D.2) 

=  ci{y)  +  jsi(y) 


with  the  aid  of  [12,  Eqs.  (3.721.2)  and  (3.721.3)].  The  sine-integral  and  eosine-integral 
may  be  eomputed  effieiently  via  rapidly  eonverging  series  representations: 


si{y)  =  -  f 

Jy 

d(y)  =  -f 

Jy 


t  2  f:f{2k-l){2k-l)\ 


cos(t)  , 
— = 

t 


C  +  ln(y)  + 


y  (-i)'y" 
4-t{2k){2ky: 


(D.3) 

(D.4) 


where  C  =  0.57721566...  is  Euler’s  eonstant.  Alternatively,  the  quantity  Ei{-jc)caxi  be 
evaluated  in  MATLAB  using  eommand  line  “cosmt(c)-A-;r/2+sinint(c))”. 
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LIST  OF  ABBREVIATIONS 


ACR 

ADC 

ADR 

AWGN 

BP 

CAN 

CDF 

CIFR 

CR 

CRN 

CS 

CSI 

DCT 

DSA 

DWT 

FDMA 

FFT 

GP 

i.i.d 

i.n.d 

FAR 

LASSO 


Adaptive  Continuous  Rate 
Analog  to  Digital  Converter 
Adaptive  Diserete  Rate 
Additive  White  Gaussian  Noise 
Basis  Pursuit 

Cooperative  Airborne  Network 

Cumulative  Distribution  Function 

Channel  Inversion  with  Fixed  Rate 

Cognitive  Radio 

Cognitive  radio  Network 

Compressive  Sampling  or  Compressed  Sensing 

Channel  Side  Information 

Discrete  Cosine  Transform 

Dynamic  Spectrum  Access 

Discrete  Wavelet  Transform 

Frequency  Division  Multiple  Access 

Fast  Fourier  Transform 

Gradient  Pursuit 

independent,  identical  and  distributive 
independent,  non-identical  and  distributive 
Least  Angle  Regression 
Least  Absolute  shrinkage  Selection  Operator 
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LB 

Lower  Bound 

LOS 

Line  of  Sight 

MIMO 

Multiple-Input-Multiple  Output 

MP 

Matching  Pursuit 

MRC 

Maximum-Ratio  Combining 

MGF 

Moment  Generating  Function 

NLOS 

Non  Line  of  Sight 

NYFR 

Nyquist  Folding  Analog-to-Information  Receiver 

OMP 

Orthogonal  Matching  Pursuit 

ORA 

Optimal  Rate  Adaptation 

OPRA 

Optimal  joint  Power  and  Rate  Adaptation 

PDF 

Probability  Density  Function 

PSD 

Power  Spectral  Density 

PU 

Primary  User 

QAM 

Quadrature  Amplitude  Modulation 

RF 

Radio  frequency 

RMSE 

Root  Mean  Square  Error 

SNR 

Signal  to  Noise  Ratio 

StOMP 

Stagewise  Orthogonal  Matching  Pursuit 

SU 

Secondary  User 

TA 

Tight  Approximation 

TDMA 

Time  Division  Multiple  Access 

TCIFR 

Eixed  Rate  with  Truncated  Channel  Inversion 
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TOMP 

UB 

WGN 


Tree  Orthogonal  Matching  Pursuit 

Upper  Bound 

White  Gaussian  Noise 
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